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Abstract 


Nonlinear  QuaJititative  Feedback  Theory  (QFT),  developed  by  Dr.  Isaac 
Horowitz,  is  used  to  design  a  flight  control  system  for  the  YF-16  aircraft.  Upon 
completing  this  stability  augmentation  system  (SAS)  additional  compensation  is 
added  to  reduce  pilot  workload  while  improving  handling  qualities. 

The  YF-16  uncertain  plant  is  simulated  with  C*  (a  blend  of  normal  accel¬ 
eration  at  pilot  station  and  pitch  rate)  as  the  controlled  output.  The  simulation 
includes  the  full  six  degree  of  freedom  nonlinear  dynamic  equations  of  motion  and 
aerodynamic  data  throughout  the  entire  subsonic  flight  envelope.  A  technique  is 
presented  which  enables  the  uncertain  nonlinear  YF-16  to  be  represented  as  a  set  of 
linear  time-invariant  plants  which  is  equivalent  to  the  nonlinear  plant  with  respect 
to  the  set  of  acceptable  outputs.  Once  this  set  of  plants  is  obtained,  a  linear  QFT 
controller  is  synthesized  yielding  fixed  compensation  which  is  extremely  insensitive 
to  varying  flight  conditions. 

Given  the  aircraft  and  resulting  flight  controller,  additional  compensation  is 
generated  which  accounts  for  the  man-in-the-loop.  The  Neal-Smith  pilot  model 
for  a  compensatory  tracking  task  is  used  to  develope  a  technique  which  allows  the 
designer  to  synthesize  compensation  which  minimizes  pilot  workload,  increases  sys¬ 
tem  bandwidth,  and  improves  handling  qualities  ratings.  The  presented  technique 
is  applicable  to  other  aircraft  and  pilot  tasks,  possessing  a  similar  pilot  model, 
given  the  desired  standards  of  performance.  This  technique  can  be  used  at  the 
earliest  stages  of  flight  control  design  phase  thus  saving  time  and  money  over  the 
current  practice. 

Simulations  in  the  time  and  frequency  domains  demonstrate  that  the  desired 
performance  is  attained.  Further  work  with  real-time  man-in-the-loop  simulations 
should  be  accomplished  to  expand  the  area  of  pilot  compensation. 


FLIGHT  CONTROLLER  DESIGN  WITH  NONLINEAR 
AERODYNAMICS,  LARGE  PARAMETER  UNCERTAINTY, 
AND  PILOT  COMPENSATION 


I.  Background 

Introduction 

The  design  of  a  flight  control  system  which  allows  the  pilot  to  maneuver  an 
aircraft  to  desired  flight  paths  is  a  major  task.  The  maneuverability  and  speeds 
required  of  modern  aircraft  place  forces  on  the  airframe  that  are  too  large  for  the 
pilot  to  control.  Also,  the  dynamic  responses  of  the  basic  airframe  are  usually 
unsatisfactory,  even  to  the  point  of  being  so  unstable  that  the  aircraft  is  totally 
uncontrollable  by  the  pilot  alone.  The  flight  controller  aids  the  pilot  in  performing 
the  tasks  which  allow  him  to  accomplish  the  mission. 

In  order  to  design  the  controller,  extensive  modeling  of  the  aircraft  dynamics 
is  required,  including  the  aerodynamic  forces.  The  parameters  in  these  models  vary 
drastically  from  one  flight  condition  to  another  and  are  not  precisely  known  for  any 
given  flight  condition.  Research  has  shown  that  the  flight  controller  design  must  be 
sufficiently  robust  to  handle  these  parameter  variations  if  the  desired  performance 
is  to  be  obtained  over  a  variety  of  flight  conditions  [19]. 

For  certain  missions  the  models  of  the  aircraft  system  are  also  highly  nonlin¬ 
ear,  with  uncertain  model  parameters  which  challenges  the  engineer  because  the 
mathematical  analysis  tools  available  are  limited.  Nonlinear  quantitative  feedback 
theory  allows  the  controls  engineer  to  take  both  of  these  aspects  into  account  at 
the  onset  of  the  design  and  guarantees  that  the  design  will  be  acceptable  for  the 
stated  problem. 


Problem  Statement 


This  thesis  uses  nonlinear  quantitative  feedback  theory  to  design  a  controller 
for  the  YF-16  allowing  the  pilot  to  control  C*,  which  is  a  handling  qualities  measure 
of  the  normal  acceleration  felt  at  the  pilot  station  blended  with  pitch  rate.  The 
controller  is  to  provide  an  acceptable  output  to  a  given  set  of  inputs  regardless  of 
the  aircraft  parameters.  In  addition  to  this  stability  augmentation  system  (SAS) 
an  outer  loop  consisting  of  compensation  to  reduce  pilot  work  load  is  designed. 
The  addition  of  this  compensation  is  to  yield  good  handling  quality  ratings  and 
increase  the  bandwidth  of  the  pilot-aircraft  system. 

Review  of  Current  Literature 

The  design  of  a  flight  controller  for  an  aircraft  is  a  problem  of  regulation 
and  control,  despite  parameter  uncertainty  and  disturbances  [8J.  The  mathemat¬ 
ical  relationships  between  the  pilot’s  inputs  and  the  output  vjiriables  are  highly 
nonlinear  as  well  as  possessing  parameters  that  are  functions  of  the  aircraft’s  flight 
conditions.  These  parameters  vary  significantly  throughout  a  flight  regime  and 
in  this  study  it  is  assumed  that  air  data  measurements  are  not  being  made  to 
determine  their  values.  These  values  would  not  be  precisely  known  even  if  such 
measurements  are  made  during  flight,  but  the  resulting  residual  uncertainty  would 
be  much  less  then  in  the  former  case.  [3).  There  is  also  the  existence  of  unknown 
external  disturbances,  such  as  wind  gusts,  whose  effect  on  the  aircraft  performance 
must  be  suitably  small  [14]. 

Current  Approach.  Presently,  there  is  no  complete  synthesis  theory  for  this 
problem.  Barfleld  [l]  stated  that  existing  flight  controllers  are  based  upon  various 
design  techniques  depending  on  the  task.  Most  of  the  controller  is  designed  around 
approximated  linear,  time-invariant  models  of  the  aircraft.  A  large  number  of 
these  models  are  then  generated  to  approximate  various  flight  conditions  which 
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the  aircraft  experiences  while  in  flight.  A  flight  controller  is  designed  for  each 
of  these  models  and  stored  in  the  on-board  memory  of  the  plane’s  computers. 
While  in  flight,  the  computer  chooses  the  controller  designed  for  the  model  that  is 
closest  to  the  current  flight  conditions  [l].  Houpis  [14]  has  indicated  that  describing 
function  theory  has  been  used  for  nonlinearities  only  to  test  for  overall  system 
stability.  The  bottom  line  is  that  current,  highly  nonlinear  aircraft  are  treated 
as  linear  systems  in  the  hope  that  the  aircraft  flight  condition  does  not  deviate 
too  far  from  the  approximations  [14].  Quinlivan  has  demonstrated  that  during 
violent  maneuvering,  the  approximations  are  violated.  The  current  solution  to 
this  problem  is  to  define  even  more  models  for  the  aircraft  during  maneuvers  in 
air-to-air  scenarios.  The  extra  models  add  to  the  number  of  controllers  to  be 
designed  to  the  point  that  different  flying  modes  must  be  selected  by  the  pilot  to 
achieve  acceptable  responses  for  different  tasks  [19'. 

Current  Techniques.  Given  the  chosen  models  to  represent  the  aircraft,  many 
techniques  exist  to  solve  the  resulting  lineaur,  time-invariant  problem.  Horowitz  and 
Shaked  concluded  that  the  superiority  of  the  transfer  function  over  the  state-  vari¬ 
able  method  explains  why  most  designs  are  accomplished  using  classical  techniques. 
Bode  plots  and  root  locus  analysis  give  insight  to  the  engineer  which  is  hidden  in 
the  state  space  representation  .9  . 

Although  the  added  insight  gained  through  the  use  of  classical  analysis  re¬ 
duces  iterations,  both  approaches  require  extensive  simulation  and  trial-and-error 
modifications.  Horowitz,  et  ai,  claimed  existing  designs,  “. .  .have  worked  because 
of  the  ingenuity  of  practical  designers  and  the  inherent  power  of  feedback;  but  a 
great  deal  of  cut  and  try  design  is  essential”  [8]. 

Uncertain  vs.  Deterministic  Models.  Horowitz,  et  al.,  suggested  that  the 
situation  is  caused  by  an  almost  total  neglect  of  quantitative  feedback  synthesis 
theory  (QFT)  in  current  designs  [8_.  Both  classical  and  modern  state-variable 
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techniques  concentrate  on  a  design  to  achieve  the  desired  output  for  a  fixed,  or 
deterministic,  model  which  represents  the  aircraft.  The  QFT  technique  provides 
a  single  design  to  give  an  acceptable  range  of  outputs  over  the  entire  range  of 
uncertainty  of  the  model  parameters  [14]. 

For  example,  consider  the  following  transfer  function  model  for  a  system: 


p(») 


K 

+  Ss  +  C 


(1.1) 


where  A,  B,  C  and  K  are  the  model’s  parameters.  Using  classical  or  modern 
techniques  requires  all  parameters  to  be  constants.  If  a  single  parameter  changes, 
a  new  design  is  required.  Using  QFT,  all  four  parameters  have  a  specified  range. 
A  design  is  generated  with  fixed  compensators  that  is  guaranteed  to  work  for 
any  combination  of  the  parameters  over  their  individual  intervals,  allowing  for  an 
infinite  number  of  possibilities  jl3]. 


The  Quantitative  Approach.  The  QFT  technique  is  summarized  in  a  techni¬ 
cal  report  [13j  for  Air  Force  Wright  Aeronautical  Labs  (AFWAL)  which  provides 
a  single  document  from  which  future  designs  can  be  referenced.  The  entire  design 
procedure  is  performed  in  the  complex  frequency  domain  to  provide  insight  as  to 
the  requirements  placed  on  the  feedback  elements  (13'.  However,  the  most  direct 
and  useful  form  of  specifications  on  a  system  are  those  on  the  time  responses. 
Krishnan  and  Cruickshanks  provided  a  simple  and  rational  means  of  translating 
time  domain  specifications  to  the  frequency  domain,  which  can  then  be  used  in 
QFT  synthesis  [17].  East  has  demonstrated  an  algorithmic  approach  to  calculate 
boundaries  on  the  loop  transmission  for  given  frequency  domain  specifications  for 
the  system.  These  specifications  include  both  the  tracking  of  a  commanded  input 
as  well  as '.attenuating  any  unwanted  disturbances  [4j. 


Houpis  suggested  a  methodological  approach  [13]  to  designing  the  controller 
once  the  above  mentioned  boundaries  are  obtained.  A  design  which  does  not 


violate  any  of  the  boundaries  is  guaranteed  to  meet  the  desired  frequency  domain 
specifications  over  the  entire  range  of  uncertainty  for  both  tracking  and  disturbance 
attenuation  [13].  Once  the  choice  of  using  quantitative  feedback  theory  has  been 
made,  there  exist  two  logical  paths  to  take  —  linear  or  nonlinear  synthesis. 


Linear  Time- Invariant  QFT.  The  simplest  application  of  quantitative  feed¬ 
back  synthesis  is  an  extension  of  the  current  flight  control  procedures;  that  is  to 
take  the  approximate  linear  models  of  the  aircraft  and  replace  the  consteint  param¬ 
eters  with  a  range  of  values.  The  resulting  designs  then  operate  at  larger  deviations 
from  the  deterministic  model. 


Although  this  technique  is  an  improvement  over  other  prominent  techniques, 
it  is  still  based  upon  the  approximate  linear  model.  However,  Hamilton  [5]  has 
shown  that  a  single  fixed  design  can  be  sufficiently  robust  to  allow  for  a  combi¬ 
nation  of  total  surface  failures,  yet  still  maintain  the  flyability  of  an  unmanned 
research  vehicle.  Simulations  of  the  design  included  nonlinearities  which  were  not 
designed  for;  however  there  was  no  guarantee  that  the  design  would  work  with  the 
nonlinear  elements  included.  Hamilton  noted  that  the  QFT  procedure  eliminates 
the  trial  and  error  portion  of  design,  but  was  still  very  tedious  due  to  the  lack  of 
computer  aided  design  programs  specifically  written  for  QFT  design  5  .  Walke  and 
Horowitz  extended  QFT  to  uncertain  MIMO  plants  whose  determinant  had  right 
half  plane  poles  and  zeros  and,  applied  it  to  the  experimental  X-29  aircraft  over 
the  entire  range  of  uncertainty  after  contracted  designers  abandoned  the  attempt 
to  independently  control  two  of  the  output  variables  [llj. 


Yaniv  has  stated  that  the  lack  of  a  generic  computer  aided  design  (CAD) 
package  has  kept  QFT  from  being  used  by  industry;  most  work  is  being  done 
by  graduate  students  [21].  Houpis  has  noted  that  AFIT  produced  ICECAP  and 
introduced  MATRIXx  packages  are  part  of  an  ongoing  effort  to  create  a  QFT  CAD 
package  but  are  not  as  yet  complete  (14|. 
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Non-Linear  QFT.  Horowitz  wsls  the  first  to  present  a  quantitative  procedure 
of  synthesis  for  nonlinear,  time-varying  systems  with  large  uncertainty.  The  work 
of  Porter  and  other  modern  theorists  at  the  time  dealt  only  with  the  qualitative 
properties  of  nonlinear  design  [10).  Since  that  time,  the  theory  has  continued  to 
expand  and  prove  itself.  Golubev  applied  the  nonlinear  QFT  technique  to  control 
a  modified  Air  Force  F-4  to  an  angle  of  attack  reaching  35  degrees,  an  extremely 
nonlinear  case  [8j. 

The  above  nonlinear  design  technique  is  based  on  conversion  of  the  nonlinear 
models  to  “equivalent  linear”  models  with  the  appropriate  uncertainty.  This  could 
lead  to  much  larger  equivalent  linear  uncertainty  than  inherently  needed.  This  is 
attributed  to  the  manner  in  which  the  “equivalent  linear"  models  were  obtained. 
Use  of  nominal  nonlinear  plant  cancellation,  involving  a  nonlinear  network  of  the 
same  complexity  as  the  plant  overcomes  this  problem  [6] .  At  a  chosen  nominal  plant 
value,  the  combination  of  the  two  behaves  like  a  linear  time-invariant  network  for  all 
inputs  and  allows  the  designer  to  greatly  reduce  over  design  [6].  The  large  “saving” 
thereby  obtained  was  vividly  demonstrated  in  a  highly  nonlinear  uncertain  2x2 
MIMO  problem  (7j. 

Barfield  points  out  the  need  for  nonlinear  uncertainty  models  in  early  flight 
control  simulations.  Current  practice  is  to  adjust  for  the  pilot  dynamics  after  the 
aircraft  is  in  the  advanced  simulation  testing  phase  where  modifications  are  expen¬ 
sive  [l].  QFT  may  prove  to  be  a  valuable  design  tool  when  working  with  nonlinear 
pilot  models  that  are  highly  uncertain.  Parameters  in  such  models  change,  not 
only  from  person  to  person,  but  from  day  to  day.  Taking  this  into  account  early 
in  the  design  may  account  for  less  need,  if  any  at  all,  of  fine  tuning  in  the  later 
stages  of  simulation. 

Summary.  The  foundation  has  been  laid  to  make  the  design  of  a  flight  control 
system  quantitative  in  nature.  The  technique  still  requires  a  large  amount  of 
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designer  interaction,  resulting  in  the  requirement  for  a  QFT  CAD  package.  The 
development  of  such  a  package  will  no  doubt  increase  the  use  of  the  technique 
throughout  the  control  industry. 

Robustness  is  inherent  in  QFT.  There  is  no  need  to  examine  system  ro¬ 
bustness  after  a  design  has  been  made,  as  is  the  case  in  much  of  modern  control 
methods. 

The  nonlinear  QFT  techniques  lend  themselves  directly  to  early  man-in-  the- 
loop  simulations  and  may  save  money  in  the  design  process. 

Finally,  it  is  obvious  that  the  modern  control  statement  that  frequency  re¬ 
sponse  methods  are  only  valid  for  linear  time- in  variant  systems,  is  false  [6,7,8,10]. 

Assumptions 

The  following  assumptions  are  made  in  the  design  problem  of  this  rhesis. 

All  parameters  contain  finite  uncertainty. 

All  commanded  inputs  and  desired  outputs  are  Laplace  transformable. 

All  aircraft  models  (plants)  have  a  unique  inverse.  This  assumption  excludes 
the  case  of  hard  saturation  of  elevator  deflection. 

The  set  of  commanded  inputs  describes  the  range  of  the  actual  inputs  to  the 
aircraft. 

The  design  is  only  in  the  longitudinal  axis,  therefore  the  only  independent 
command  input  is  the  elevator. 

It  should  be  noted  that  assumptions  3  and  5  are  added  here  only  to  reduce  the 
scope  of  the  thesis  and,  in  theory,  can  be  releixed.  QFT  has  been  extended  to 
handle  the  hard  saturation  problem,  but  the  usual  QFT  can  be  applied  only  if  the 
inverse  range  is  finite.  For  example  given  an  output  y(f),  there  is  only  a  finite 
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number  of  inputs  i,(t)  which  can  generate  this  y{t)  [8j.  In  actuality,  uncertainty 
can  be  semi-infinite,  i.e.  a  gain  factor  from  1  to  oo  but  not  from  — oo  to  +oo  [8]. 

Scope 

The  foundation  of  this  thesis  is  laid  in  nonlinear  quantitative  feedback  theory. 
This  thesis  applies  the  technique  to  a  single  input,  single  output  (SISO)  nonlinear 
system  with  large  parameter  uncertainty.  There  is  an  inner  stability  augmenta¬ 
tion  system  (SAS)  loop  involving  a  nonlinear  airframe  model  and,  an  outer  loop 
containing  the  pilot,  where  the  model  contains  a  pure  time  delay. 

The  full,  nonlinear  six  degrees  of  freedom  equations  of  motion  are  used  to 
generate  the  plants,  although  the  design  is  done  only  for  the  longitudinal  axis. 

The  only  commanded  input  to  the  system  is  a  symmetric  elevator  deflection 
and  the  commanded  output  is  C*:  normal  acceleration  at  the  pilot  station  blended 
with  pitch  rate. 

The  completed  design  is  simulated  as  a  linear  system  using  MATREXx  and 
on  a  Fortran  simulator  using  the  six  degrees  of  freedom,  nonlinear  equations  of 
motion.  The  validation  of  the  design  is  to  demonstrate  that  all  of  the  output 
responses  are  within  acceptable  limits. 

Approach 

The  first  step  in  the  design  procedure  is  to  obtain  linear  time- invariant  models 
which  are  rigorously  equivalent  to  the  nonlinear  plant  with  respect  to  the  defined 
set  of  desired  outputs.  These  equivalent  plants  are  generated  from  time  histories  of 
the  input  and  output  from  simulator  data  supplied  by  AFWAL.  This  data  is  used 
in  a  program,  developed  by  Golubev  [8]  and  provided  by  Dr.  Horowitz,  to  generate 
linear  time-invariant  (LTI)  plants  that  are  equivalent  to  the  nonlinear  plant  in  that 
both  models  give  the  same  output  for  the  given  input.  This  procedure  is  repeated 
for  the  set  of  inputs  and  outputs  that  are  to  exist  in  the  operating  region  of  the 
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system.  The  result  is  a  set  of  LTI  equivalent  plants  which  are  used  in  the  design 
process. 

Next,  the  boundaries  on  acceptable  outputs  provided  by  AFWAL  are  trans¬ 
lated  from  the  time  domain  to  the  frequency  domain.  At  this  point,  a  decision  can 
be  made  whether  the  specifications  are  reasonable  and  attainable  by  any  design 
method.  QFT  or  other.  The  generation  of  bounds  on  a  nominal  loop  in  the  fre¬ 
quency  domain  are  obtained,  which  must  be  satisfied  in  order  that  the  output  is  in 
the  acceptable  set  over  the  entire  set  P .  When  the  design  is  completed  without  vio¬ 
lating  any  boundary,  the  controller  is  guaranteed  to  work,  not  only  for  the  nominal 
plant,  but  also  for  the  entire  nonlinear  plant  set  over  its  range  of  uncertainty. 

The  controller  design  is  linearly  simulated  on  MATRDCx  as  a  first  validation 
and  a  design  tool.  The  final  validation  is  performed  on  an  existing  nonlinear  YF-16 
simulator  using  the  full  six  degree  of  freedom  equations. 


Presentation 


This  thesis  has  seven  chapters.  Chapter  II  develops  the  nonlinear  QFT  the¬ 
ory.  Chapter  III  presents  the  generation  of  the  linear  time- in  variant  plants.  Chap¬ 
ter  IV  explains  the  development  of  the  SAS  compensation  of  the  inner  loop.  Chap¬ 
ter  V  presents  the  MATRIXx  linear  and  Fortran  nonlinear  simulations  of  the  inner 
loop  excluding  pilot  compensation.  Chapter  VI  explains  the  development  of  the 
additional  single  degree  of  freedom  pilot  compensation,  as  well  tis  its  simulation, 
while  Chapter  VII  contains  the  conclusions  and  recommendations  for  future  study. 
Appendix  A  includes  the  LTI  set  P  derived  from  Golubev’s  program  along  with 
plots  validating  these  transfer  functions.  Appendix  B  presents  the  plant  templates 
of  interest  and  frequency  domain  tracking  bounds.  Appendix  C  contains  simula¬ 
tions  of  the  inner  loop  SAS  along  with  the  data  used  to  obtain  the  prefilter,  while 
Appendix  D  derives  the  fourth  order  Fade  approximation  of  a  pure  time  delay  used 
in  the  Fortran  simulations. 
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II.  Nonlinear  QFT  Theory 


Introduction 

This  chapter  presents  the  theory  behind  controller  design  for  a  nonlinear 
uncertain  plant.  A  brief  overview  of  linear  SISO  QFT  design  is  presented.  The 
reader  is  referred  to  a  technical  report  [13]  for  a  comprehensive  presentation  of  the 
design  process. 

Linear  SISO  Case 

In  the  design  of  a  controller  for  a  linear  model  the  technique  begins  with  a 
set  of  linear  differential  equations  whose  parameters  are  uncertain,  giving  a  set  of 
plant  transfer  functions  P  ~  {P}.  Every  possible  combination  of  parameters  gives 
a  different  plant  contained  in  the  set  P. 

It  is  desired  that  the  final  closed  loop  tracking  response  be  within  given 
or  desired  acceptable  bounds,  for  any  possible  combination  of  parameters  and 
therefore  any  P  €  F .  If  the  acceptable  bounds  are  given  in  the  time  domain  they 
can  be  converted  to  equivalent  frequency  domain  bounds  !l3^  Once  the  closed 
loop  frequency  response  bounds  are  obtained,  the  compensator  G]^)  of  Figure  2.1 
is  used  to  insure  that  the  variation  in  the  closed  loop  response  is  less  than  or  equal 
to  that  of  the  bounds.  However,  this  is  insufficient  to  guarantee  that  the  closed  loop 
tracking  response  for  every  P  £  P  lies  within  the  tracking  bounds.  The  prefilter 
P(s)  of  Figure  2.1  is  used  in  this  final  task. 

The  designs  of  the  compensators  G[s)  and  F[$)  are  based  upon  boundaries 
that  are  derived  for  a  set  of  frequencies.  Similar  bounds  on  G{s)  can  be  obtained  for 
the  response  of  the  system  to  external  disturbances.  All  of  the  bounds  mentioned 
above  can  be  generated  manually  on  the  Nichols  chart  [13]  or  the  same  concepts 
can  be  programed  on  a  computer  for  automation  as  described  in  Chapter  IV'. 
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Figure  2.1.  QFT  Compensated  Block  Diagram 


Nonlinear  SISO  Case 

The  problem  with  the  nonlinear  case  is  at  the  onset  of  the  design.  The 
differential  equations  are  nonlinear  and  therefore  transfer  functions  may  not  be 
used  in  the  usual  sense.  There  exists  a  set  of  nonlinear  plants  Ti'  due  to  every  ^ 

possible  combination  of  parameter  values.  It  is  still  desired  that  all  system  outputs 
due  to  realistic  inputs  lie  within  an  acceptable  set  A. 

It  is  necessary  to  find  a  set  of  linear  time-inveuiiint  plants  P  that  is  precisely  ^ 

equivalent  to  the  nonlinear  set  "V)  with  respect  to  the  acceptable  output  set  A.  A 
rigorously  proven  technique  for  this  problem,  using  Schauder’s  fixed  point  theo¬ 
rem,  shows  that  a  design  meeting  specifications  for  the  set  P  will  also  satisfy  the  | 

nonlinear  problem  with  W  ilOj. 

Direct  Analytical  Model.  For  the  nonlinear  plant  of  Figure  2.2  having  the 
input  x(t),  and  the  output  y{t)  is  of  the  form; 
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y{t)  =  w{x{t)) 


(2.1) 


Given  any  Ci  €.  A  and  any  plant  €  W,  then  the  problem  becomes  one  of 
determining  the  input  to  which  gives  the  output  Cj  by  solving  Eq  (2.1)  for  i 
given  y.  The  resulting  i  is  denoted  by  xj.  It  is  necessary  to  find  a  linear  time- 
invariant  plant  P/  which  is  equivalent  to  for  this  one  case.  That  is,  the  output  of 
P/  is  a,(f)  when  its  input  is  x((<).  When  a, (t)  and  x^(t)  are  Laplace  transformable, 
a  difficult  condition  to  violate,  one  can  write  [8!: 


Laijt) 


(2.2) 


The  linear  time-invariant  P/(s)  is  precisely  equivalent  to  the  nonlinear  only 
when  the  input  is  .  In  order  to  describe  the  nonlinear  set  W  with  a  linear  set  P , 
a  set  of  outputs  and  corresponding  inputs  is  required.  This  can  be  accomplished 
by  repeating  the  above  technique  over  all  the  w  ^  W ,  for  the  same  output  o,(t). 
Then  repeating  the  above  steps  for  all  of  the  desired  outputs  a  A,  generates  a 
set  of  {/V}  =  P  [8]. 
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Non-Direct  Analytical  Model.  Another  technique,  based  in  theory  on  the 
above  procedure,  is  very  useful  for  nonlinear  plants  that  have  extremely  compli¬ 
cated  analytic  models.  This  technique  requires  either  the  physical  availability  of 
the  plant  or  a  nonlinear  simulation  of  the  plant.  In  this  thesis  a  simulation  of 
the  YF-16  aircraft  using  the  nonlinear  equations  of  motion,  complete  with  a  large 
data  base  of  aerodynamic  data,  is  used.  A  stabilizing  controller  is  first  obtained, 
by  “intelligent”  cut  and  try,  around  the  plant  with  the  only  goal  being  that  the 
acceptable  range  of  outputs  can  be  obtained;  no  attempt  is  made  at  a  reasonable 
design.  The  simulator  is  then  run,  with  the  operator  supplying  the  “input”  as 
needed  to  obtain  one  acceptable  output.  A  record  is  kept  of  the  input  signal  to  the 
plant  and  the  plant’s  output.  This  can  be  done  by  writing  the  two  time  histories 
to  data  files.  The  simulation  is  repeated  until  a  set  of  acceptable  outputs  is  ob¬ 
tained  that  fills  the  envelope  of  acceptable  responses.  The  set  of  plant  output  and 
corresponding  input  time  history  points  is  then  used  to  generate  the  set  of  linear 
equivalent  plants  P. 

Summary 

This  chapter  presents  an  overview  of  the  continuous  SISO  QFT  design  tech¬ 
nique.  with  extension  to  the  nonlinear  case.  Two  techniques  for  obtaining  the 
linear  time-invariant  equivalent  plants  are  presented,  either  of  which  may  be  used 
to  accomplish  the  task. 
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III.  Derivation  of  the  LTI  Equivalent  Plants 


Introduction 

As  noted  in  the  previous  chapter  this  thesis  does  not  solve  the  analytical 
model  backwards  to  obtain  the  linear  time-invariant  equivalent  plants.  Instead  a 
simulator  provided  by  Tom  Cord  of  AFWAL  is  used.  The  simulator  was  originally 
developed  for  studying  the  high  angle  of  attack  and  the  spin  characteristics  of  the 
YF-16. 


The  Simulator 

The  following  modifications  are  made  to  the  original  simulator  in  order  to 
aid  in  obtaining  the  desired  data: 

Longitudinal  stick  force  shaping  compensation  is  removed. 

The  commanded  output  C*  is  generated  as  an  output  variable. 

Various  output  files  for  transfer  function  generation  and  MATRIXx  graphics 
are  added. 

Normal  acceleration  feedback  is  replaced  with  C’  feedback. 

Prefiltering  of  input  signals  is  added. 

The  simulator  is  a  Fortran  implementation  of  the  actual  control  laws  as  well 
as  the  six  degree  of  freedom  nonlinear  equations  of  motion  for  the  aircraft.  The 
aerodynamic  derivatives  for  the  equations  of  motion  are  supplied  from  data  tables 
covering  the  entire  fiight  envelope,  including  high  smgle  of  attack.  During  the 
simulation,  the  equations  of  motion  are  solved  using  a  Runge-Kutta  integration 
scheme  with  the  states  updated  every  5  msec.  Reduction  of  the  update  interval 
shows  no  improvement  in  the  accuracy  obtained  as  compwed  to  actual  fiight  test 
data. 
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Each  time  the  simulation  is  started,  a  prompt  for  load  factor  in  g’s  aind  angle 
of  attack  in  degrees  is  used  to  trim  the  aircraft  to  its  nominal  flight  condition.  This 
thesis  has  the  nominal  condition  of  straight  and  level  flight  at  0.9M  and  20,000  ft 
in  the  derivation  of  the  equivalent  plants. 

The  input  to  the  plant  is  considered  to  be  the  sum  of  the  actuating  signals 
to  the  two  elevators  of  Figure  3.1,  which  are  always  symmetric  since  the  elevator 
is  the  only  commanded  input.  The  output  is  C*  in  g’s  which  for  the  YF-16  is 

C‘  =  +  12.4g  (3.1) 

where  An^,  is  the  normal  acceleration  at  the  pilot  station  in  g’s,  and  q  is  the  pitch 
rate  in  rad/sec.  By  specifying  the  above  input  and  output  the  plant  is  defined  as  the 
airframe  with  actuators,  including  the  position  and  rate  limiters  of  the  actuators. 
Also  included  are  the  dynamics  of  the  aircraft’s  leading  edge  flaps  of  Figure  3.1. 
Since  this  thesis  is  restricted  to  a  single  input,  the  two  available  choices  are  to  either 
disable  the  flaps  entirely  in  the  trailing  position,  or  make  them  a  dependent  input. 
The  latter  is  chosen  by  commanding  the  leading  edge  flaps  with  an  angle  of  attack 
plus  pitch  rate  schedule.  This  not  only  gives  an  input  totally  dependent  upon  the 
commanded  elevator  input,  but  also  has  a  stabilizing  effect  on  the  aircraft. 

The  Equivalent  Plants 

In  order  to  adequately  represent  with  the  set  P  it  is  necessary  to  know 
the  range  of  acceptable  outputs  a  priori  and  drive  the  simulator  to  outputs  that 
adequately  represent  that  envelope.  In  Figure  3.2  are  the  proposed  normalized 
bounds  on  a  C*  response  to  a  step  input  used  in  this  study.  The  lower  boundary 
is  modified  from  that  proposed  by  Tobie  [20].  The  bottom  of  this  envelope  does 
not  allow  a  higher  order  response  and  specifies  a  minimum  greater  than  zero  for 
the  first  0.2  seconds.  This  modification  makes  the  lower  bound  more  severe  then 
the  original.  The  purpose  behind  the  modification  is  to  provide  a  smoother  lower 
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Figure  3.1.  YF-16  and  Control  Surfaces 


bound  on  the  closed  loop  response  in  the  frequency  domain,  as  well  as  to  disallow 
non-minimum  phase  responses,  with  no  sacrifice  to  the  time  response. 

The  objective  of  the  simulations  is  to  drive  a  set  of  outputs  that  fill  the 
envelope  of  Figure  3.2.  This  is  accomplished  with  8  outputs  for  a  Ig  envelope 
and  7  cases  for  a  2g  envelope,  all  15  normalized  cases  are  shown  in  Figure  3.3 
superimposed  upon  the  normalized  bounds  and  are  judged  to  populate  them  very 
well. 

The  next  objective  is  to  find  15  linear  time- invariant  trzinsfer  functions,  one 
for  each  input/output  pair,  given  the  time  histories  obtained  from  the  simulator. 
Golubev’s  Fortran  subroutine  provided  by  Dr.  Isaac  Horowitz  is  used  to  help  ac¬ 
complish  this  task.  By  providing  a  vector  of  the  plant’s  input,  output,  and  time 
of  occurrence  the  subroutine  utilizes  the  Numerical  Algorithms  Group  (NAG)  For¬ 
tran  Library  and  Grout's  factorization  method  to  obtain  a  linear  transfer  function 
equivalent  to  the  nonlinear  plant  for  that  one  specific  input.  The  subroutine  re- 


3-3 


TIME  IN  SECONDS 


Figure  3.2.  Normalized  Bounds  on  the  Output  C* 
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Figure  3.3.  Simulator  Outputs  used  to  Generate  Plants 
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quires  the  user  to  specify  the  order  of  both  the  numerator  and  denominator;  a 
small  amount  of  experience  in  analyzing  the  time  responses  of  transfer  functions 
allows  the  user  to  eliminate  most  of  the  trial  and  error  in  finding  a  suitable  fit.  The 
decision  is  made  to  fit  the  time  histories  with  transfer  functions  that  all  have  the 
same  number  of  excess  poles  over  zeros,  while  placing  no  constraint  on  the  order 
of  the  system  model.  Thus  causing  the  plant  templates  to  have  no  uncertainty  in 
phase  angle  at  w  =  oo.  Analyzing  the  initial  10  msec  of  the  response  gives  good 
insight  into  selection  of  the  actual  number  of  excess  poles.  An  excess  of  two  poles 
is  chosen. 

In  addition  to  validating  the  program  with  the  input/output  time  histories  of 
known  transfer  functions,  known  test  cases  are  run  on  the  simulator.  The  simulator 
is  given  a  step  input  with  the  leading  edge  flaps  frozen  and  the  output  recorded. 
This  data  is  then  fitted  with  a  transfer  function  that  matches  that  of  a  linearized 
YF-16,  without  leading  edge  flaps,  at  the  same  nominal  flight  condition.  The 
final  trzmsfer  function  fits  are  validated  by  reapplying  the  original  input  to  the 
proposed  transfer  function  and  comparing  this  output  with  the  original  output. 
The  two  outputs  are  then  plotted  along  with  their  error.  Excellent  results  are 
obtained  as  shown  in  Figure  3,4.  The  results  of  all  15  cases  along  with  the  transfer 
functions  obtained  are  contained  in  Appendix  A. 

Summary 

This  chapter  presents  the  derivation  of  the  set  of  plant  linear  time-invariant 
equivalent  transfer  functions  used  to  represent  the  nonlinear  YF-16.  The  simulator 
is  described  in  detail,  as  well  as  the  modifications  made  to  it.  The  tracking  response 
bounds  on  which  the  final  design  is  evaluated,  are  presented  for  the  first  time.  The 
use  of  the  response  bounds  in  operating  the  simulator  is  explained  and  examples 
presented.  Converting  time  histories  into  transfer  functions,  along  with  a  validation 
technique  are  both  presented. 
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IV.  Design  of  the  Inner  Loop  Controller 


Introduction 

This  chapter  discusses  the  design  of  the  compensators  F{s)  and  G(s)  shown 
in  Figure  2.1.  The  bounds  on  the  nominal  loop  transmission  La(s)  are  modified 
slightly  from  those  obtained  in  the  usual  QFT  problem,  due  to  different  responses 
obtained  from  Ig  and  2g  commands.  After  the  modified  bounds  are  established, 
the  nominal  loop  is  shaped  on  the  Nichols  chart.  Once  the  closed  loop  is  completed, 
bounds  on  the  prefilter  F(s)  are  derived  and  the  compensator  designed. 

Bounds  on  Lo{s) 

The  first  step  in  generating  the  bounds  on  the  nominal  loop  transmission 
is  converting  the  time  domain  bounds  of  Figure  3.2  into  bounds  on  the  closed 
loop  response  in  the  frequency  domain  as  shown  in  Figure  4.1.  This  is  usually 
performed  using  standard  figures  of  merit  along  with  some  trial  and  error;  however, 
the  program  used  to  generate  the  equivalent  plants  can  also  be  used. 

In  obtaining  the  equivalent  plants,  it  is  noted  that  the  2g  responses  alone  do 
not  populate  an  area  of  the  envelope  very  well  due  to  inertia  of  the  aircraft  and  use 
of  the  same  bounds  for  Ig  and  2g  commands.  To  make  the  amount  of  uncertainty 
more  realistic  an  additional,  more  severe,  upper  bound  is  placed  on  the  normalized 
2g  responses  (see  Appendix  B).  Since  for  all  cases  the  output  is  a  minimum  phase 
response,  only  the  magnitudes  of  the  frequency  domain  bounds  are  required  to 
generate  the  bounds  on  Z-o(s). 

The  transfer  functions  used  to  model  the  Ig  upper  bound,  2g  upper  bound, 
and  common  lower  bound  are  respectively: 


15.8965(5  +  0.9022) 
(s^  2.6342)(5  ^  5.439) 


(4.1) 


Figure  4.1.  Closed  Loop  Tracking  Bounds  in  the  Frequency  Domain 


7,4435(a  +  1.5262) 

(5  +  2.112  ±  >2.5852) 

19.2109 

(.s  +  1.5942)  (s  +  12.088) 

For  every  frequency  there  is  a  maximum  allowable  variation  in  the  magnitude 
of  the  response  due  to  the  set  of  plants  P.  The  change  in  the  calculation  of  the 
bounds  on  Lo(ju})  mentioned  earlier  is  handled  in  the  manner  described  next. 

Modified  Tracking  Bounds.  A  plant  template,  consisting  of  the  magnitude 
and  phase  angle  of  each  of  the  15  plants,  at  a  frequency  u  =  a;,-,  is  generated 
in  the  usual  manner  [l3j.  However  on  this  combined  template,  two  subtemplates 
should  also  be  distinctly  marked,  along  with  a  single  nominal  plant’s  magnitude  and 
phase.  One  subtemplate  consist  of  only  the  points  due  to  the  Ig  response  (which 
contains  the  nominal  plant),  while  the  other  is  for  the  2g  response.  Examples  of  two 
combined  templates,  each  representing  a  different  magnitude  but  the  same  phase 
angle  of  the  nominal  plant,  are  shown  in  Figure  4.2.  The  Ig  subtemplate  consists 
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Figure  4.2.  Example  of  Modified  Tracking  Bounds  at  '*;  =  1  rad/sec 


of  the  eight  plant  points  between  A  (A')  and  B  [B')  while  the  2g  subtemplate 
consists  of  the  seven  plant  points  between  C  {€')  and  D  {D'). 

The  combined  Ig  and  2g  template  [ABCD]  is  then  placed  on  the  Nichols 
chart  to  locate  a  point  on  the  boundary  of  Lo{jui),  using  the  allow'able  variation 
in  magnitude  from  the  Ig  bounds  in  accordance  with  the  technique  for  obtaining 
tracking  bounds  given  by  linear  continuous  QFT.  For  the  remainder  of  this  thesis 
unmodified  QFT  techniques  are  referred  to  as  the  usual  case.  The  usual  bound 
at  —100°  for  1  rad/sec  is  given  by  the  position  of  the  template  ABCD  in  Fig¬ 
ure  4.2  where  the  nominal  point  is  indicated  by  x  .  The  maximum  closed  loop 
magnitude  of  1.399  dB  occurs  at  point  C  and  the  minimum  of  -2.864  dB  occurs 
at  point  B,  giving  a  variation  of  4.263  dB  for  the  combined  Ig  and  2g  template. 
The  allowable  variation  at  1  rad/sec  for  the  Ig  bound  is  4.264.  The  allowable 
variations  can  be  obtained  from  Figure  4.1  or  Tables  B.l  and  B.2.  Thus,  the  nom¬ 
inal  point  of  the  template  would  be  marked  as  a  point  on  the  bound  of  Lo(jl); 
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Figure  4.3.  Calculation  of  Bounds  on  F(jl)  for  the  ABCD  Template 

however,  instead  of  marking  the  nominal  plant’s  position,  one  more  check  must  be 
performed.  Note  the  maximum  and  minimum  closed  loop  magnitudes  due  to  the 
Ig  subtemplate  only,  in  Figure  4.2  these  are  -0.876  db  at  point  A  and  -2.864  dB 
at  point  B.  Compare  these  two  extreme  values  with  the  maximum  and  minimum 
allowed  b}'  the  Ig  tracking  response  bound  of  2.766  dB  and  -1.497  dB  from  Fig¬ 
ure  4.1.  Next  calculate  the  range  of  additional  magnitude  required  to  meet  the 
bounds  at  w  =  1.  The  maximum  gain  that  cam  be  added  by  F{jl)  to  lie  within 
the  bounds  of  Figure  4.1  is  2.766  —  (—0.876)  =  3.642  dB  while  the  minimum  is 

—  1.497  —  (—2.864)  =  1.367  dB  sis  illustrated  in  Figure  4.3.  Thus  the  range  of  al¬ 
lowable  gain  increase  is  from  1.367  dB  to  3.642  dB.  This  range  of  magnitudes  is 
the  bound  on  the  prefilter  F(jl)  if  the  nominal  loop  is  actually  at  1.47  dB  and 

—  100°  at  w  =  1  rad/sec,  i.e.  Lo{jl)  =  1.47dBZ  -  100°. 

Next  repeat  the  same  calculation  on  the  bounds  of  F{jl)  except  this  time 
use  the  2g  subtemplate  and  2g  closed  loop  tracking  bounds.  For  the  2g  subtemplate 
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the  maximum  closed  loop  magnitude  is  1.399  dB  at  point  C  while  the  minimum  is 
-0.404  dB  at  point  D.  Figure  4.1  (also  see  Table  B.2)  gives  the  maximum  allowable 
closed  loop  magnitude  1.844  dB  and  the  minimum  of  -1.497  dB.  Calculating  the 
acceptable  bounds  on  F{jl)  gives  a  maximum  of  1.844  —  1.399  =  0.445  dB  and 
a  minimum  of  —1.497  —  (—0.404)  =  —1.093  dB.  Compare  the  two  sets  of  bounds 
on  F{jl)  thus  obtained.  The  goal  of  this  check  is  to  have  one  prefilter  F{juj)  to 
meet  system  specifications.  Therefore,  if  there  is  some  overlap  between  the  two 
sets  of  bounds  on  F(ju,’i)  there  exists  a  single  F{jLj)  which  meets  both  bounds, 
and  the  boundary  location  is  a  valid  boundary  point  of  Lo{u;i)-  If  there  is  no 
overlap  between  the  two  sets  of  bounds  on  F{ju!,)  then  the  open  loop  gain  must 
be  increased  at  this  phase  value,  moving  the  combined  template  up  the  Nichols 
chart  until  an  acceptable  point  is  located.  For  the  template  ABCD  in  Figures  4.2 
and  4.3  the  ranges  of  3.642  dB  to  1.367  dB  and  0.445  dB  to  -1.093  dB  have  no 
overlap.  This  is  graphically  depicted  in  Figure  4.3  where  for  the  Ig  subtemplate 
the  minimum  increase  in  gain  of  1.367  dB  is  greater  then  the  maximum  allowable 
of  0.445  dB  for  the  2g  subtemplate.  Therefore,  the  gain  of  the  plant  at  -100®  must 
be  increased  raising  the  combined  template. 

If  the  nominal  point  of  the  template  is  raised  to  3.82  dB  the  template 
A'B'C'D'  of  Figure  4.2  is  obtained.  Since  the  allowable  variation  is  met  at  the  previ¬ 
ous  lower  gain  no  additional  check  is  required.  Points  A'  and  B'  have  the  closed  loop 
magnitudes  of  -0.134  dB  and  -2.09  dB  respectively  giving  an  acceptable  Ig  range 
on  |F(il)l  of  2.766  -  (-0.134)  =  2.9  dB  to  -1.497  -  (-2.09)  =  0.593  dB.  Points 
C  and  D'  have  the  magnitudes  1.248  dB  and  -0.241  dB  respectively  giving  an  ac¬ 
ceptable  2g  range  on  IF(jl)  of  1.844  -  1.248  =  0.596  dB  to  -1.497  -  (-0.241)  = 
—1.256  dB.  There  now  exists  an  overlap  between  the  ranges  on  |F(jT)l,  thus 
3.82  dB  at  —100°  is  a  valid  boundary  point  of  Lo{jl).  This  process  is  repeated 
at  various  values  of  open  loop  phase  angle  for  u;  =  1  rad/sec  and  then  for  various 
values  of  Wj. 
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Finding  the  closed  loop  magnitude  M,  for  the  form  of  Figure  2.1  without 
the  prefilter  F{s),  when  given  the  open  loop  magnitude  m  and  phase  angle  6,  i.e. 
L{ju>)  =  ml9,  can  be  easily  automated  by  solving  for  the  magnitude  of  the  closed 
loop  transfer  function: 


lM(ja;)|  = 


1  +  L{juj) 


When  both  M  and  m  are  in  arithmetic  units  not  decibels,  then 


V^l  +  +  2m  cos  6 

This  method  is  used  to  obtain  the  magnitudes  of  the  above  example  and  through 
the  remainder  of  this  thesis. 


Stability  Bounds.  The  modified  tracking  bounds  are  obtained  for  templates 
up  to  a  frequency  for  which  the  allowable  closed  loop  variation  is  sufficiently  greater 
then  the  variation  in  the  templates’s  magnitude.  In  this  thesis  a  frequency  of 
Wh  =  10  rad/sec  is  sufficient.  Above  this  frequency,  only  stability  bounds  exist 
representing  the  maximum  allowable  closed  loop  magnitude  [13].  This  thesis  uses 
a  stability  or  M  contour  of  3  dB.  The  boundaries  on  the  nominal  loop  transmission, 
Lo{ju)  are  shown  on  the  Nichols  chart  in  Figure  4.4. 

The  Shaping  of  Lo{s] 

Next  the  nominal  loop  Lo{ju})  is  shaped  by  choosing  an  appropriate  G{s) 
such  that  Lo{ju>i)  remains  outside  the  boundary  B(u;,)  for  each  w,  =  w.  Since 
Lol-s)  =  G(s)/*o('S)-  This  is  accomplished  by  choosing  the  poles,  zeros,  and  gain 
of  G{s).  The  optimal  loop  transmission  Xopt  lies  on  B(uj)  at  all  u>  and  not  only 
exists  but  is  unique  [81.  The  design  of  a  practical  I<o(s)  is  somewhat  of  an  art. 
The  greater  the  number  of  poles  and  zeros  of  G(s)  that  are  chosen,  the  closer  to 
the  optimum  Lo(s)  is  realized,  thus  there  e.xists  a  tradeoff  between  complexity  and 
bandwidth  l8l. 
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Open  Loop  Phase  (deg) 


Template  Magnitude  Variation 
Greater  Than  Allowed  by  Bounds 

Template  Magnitude  Variation 
Less  Than  Allowed  by  Bounds 

1  Frequency  Ranges 

0  <  w  < 

(jJl  <  UJ  <  Ufh 

u>h  <  u  <  oo 

j  on  or  above 

on  or  above 

j  and 

and 

stability  only 

'  stability 

on  or  below 

Table  4.1.  Frequency  Ranges  of  Various  Boundaries 


Interpretation  of  the  Boundaries.  The  reader  is  cautioned  on  the  interpre¬ 
tation  of  “remaining  outside  the  boundaries” .  Most  of  the  literature  breaks  down 
the  bounds,  B{s),  into  two  subsets  as  is  done  here.  However,  one  may  be  led  to 
believe  that  the  loop  transmission  must  be  outside  the  stability  bounds  but  “on  or 
above”  the  tracking  bounds.  This  is  not  necessarily  the  case.  There  exits  a  range 
of  frequencies  where  Lo[j'^)  may  have  to  be  below  a  certain  boundary  to  meet 
tracking  bounds.  The  boundary  B{8)  in  Figure  4.4  shows  this  exception.  In  cases 
where  lLo{j8)  >  -125°,  the  magnitude  of  Lo{j8)  must  be  greater  than  the  point 
on  the  B{8)  boundary  at  that  phase.  This  is  the  normal  case,  Lo{j8)  must  be  on 
or  above  B{8)  on  the  Nichols  chart.  How'ever,  for  lLo{j8)  <  -125°  the  magnitude 
must  be  less  than  the  point  on  the  boundary  at  that  phase  to  insure  that  the 
closed  loop  variation  in  magnitude  is  less  than  that  allowed  at  u.'  =  8  rad/sec  in 
Figure  4.1.  Therefore  there  exists  a  region  where  Lo[jui)  must  be  “on  or  below” 
B{u}).  It  is  emphasized  that  the  entire  boundary  B(8)  is  a  tracking  bound,  not 
partly  tracking  bound  and  partly  stability  bound  as  in  the  czise  for  B(2.5). 

The  range  of  frequencies  of  interest,  where  one  must  check  for  the  on  or 
below  condition,  begins  at  the  frequency  w  =  Wj  where  the  magnitude  variation 
of  the  template  first  becomes  less  than  the  variation  allowed  by  the  closed  loop 
frequency  bounds  of  Figure  4.1.  The  range  ends  at  =  u.>.  The  definition  of  this 
range  results  in  a  quantitative  method  of  obtaining  ujh  as  the  highest  frequency 
for  which  there  exists  an  on  or  below  region  of  the  tracking  boundary.  Therefore, 
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for  frequencies  above  Wh  only  stability  bounds  exist.  These  frequency  ranges  are 
summarized  in  Table  4.1.  A  quick  method  to  check  if  a  frequency  has  an  on  or 
below  boundary  is  to  place  the  template  at  the  bottom  of  the  M  contour  on  the 
Nichols  chart,  as  if  determining  stability  bounds,  and  checking  the  actual  variation 
in  closed  loop  magnitude.  Since  the  closed  loop  contours  are  the  closest  together 
near  the  0  dB,  —180°  point,  then  this  region  on  the  chart  is  where  the  largest 
variation  in  the  closed  loop  magnitude  exists  while  meeting  the  stability  bound.  If 
the  closed  loop  variation  is  satisfied,  the  stability  bound  dominates.  If  the  closed 
loop  variation  is  not  met,  the  template  on  the  Nichols  chart  must  be  lowered  until 
the  variation  is  satisfied  and  the  resulting  boundary  is  an  “on  or  below”  boundary. 
For  this  thesis  the  region  to  check  is  approximately  3  <  w  <  10  rad/sec. 

Proposed  G(s).  The  first  design  of  G(s)  is  accomplished  using  the  bounds 
on  Loijui)  of  Figure  4.4  on  page  4-7  with  the  exception  of  B(18),  and  is  shown 
superimposed  upon  these  bounds.  A  prefilter  F(s)  is  then  designed  in  the  usual 
manner  (13).  This  design  is  then  examined  using  linear  simulation  on  MATRDCx 
by  obtaining  the  responses  of  each  of  the  15  LTI  plants  in  P  using  the  fixed  G{s) 
and  F{s).  All  responses  are  within  the  specified  boundaries  and  all  except  one 
are  judged  to  be  satisfactory.  The  response  which  uses  plant  15  has  a  lightly 
damped  sinusoidal  component  of  approximately  16  rad/ sec  superimposed  upon  an 
acceptable  first  order  response.  Plots  of  the  closed  loop  frequency  responses  give 
immediate  insight  into  the  problem  as  can  be  seen  in  Figures  4.5  and  4.6. 

The  decision  is  made  to  proceed  to  the  nonlinear  simulation  based  upon  the 
following  reasoning.  The  commanded  inputs  to  the  simulator  that  produced  the 
time  histories  to  form  plant  15  are  extremely  erratic,  worse  than  any  other  case.  It 
is  felt  that  this  plant  represents  a  very  unlikely  operating  condition  for  the  aircraft 
and  is  therefore  ignored  for  the  time  being.  The  error  in  this  reasoning  is  that, 
the  signal  to  be  analyzed  in  making  the  judgement  of  erratic  response  is  the  input 
signal  to  the  aircraft  actuators  which,  in  fact,  is  rather  well  behaved.  The  nonlinear 
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Figure  4.5.  Response  of  Plant  #  15  for  Initial  Design 


Figure  4.6.  Closed  Loop  Frequency  Response  of  Plant  #15 


simulation  points  this  out  by  exciting  a  limiting  cycle,  while  attempting  to  trim,  at 
»  16  rad/sec.  This  agrees  with  the  resonant  peaks  seen  in  the  linearly  simulated 
frequency  response  of  plant  15  in  Figure  4.6.  A  closer  examination  of  the  bounds 
for  15  <  0/  <  20  rad/sec  from  which  the  boundary  B(18)  in  Figure  4.4  is  obtained, 
makes  very  clear  that  this  boundary  is  violated.  This  first  design  points  out  the 
drawbacks  of  representing  the  infinite  number  of  boundaries  with  a  finite  set,  but 
exemplifies  the  insight  possible  from  working  in  the  frequency  domain. 

Final  Design  of  G(s).  Using  the  full  set  of  bounds,  JB(5),  on  Lo(s)  of  Fig¬ 
ure  4.4  a  new  nominal  loop  transmission  is  designed.  The  emphasis  again  is  on 
minimizing  the  bandwidth  of  G(s)  while  having  Loi-s)  remain  outside  of  the  bound¬ 
aries,  and  on  accomplishing  this  with  a  practical  G(s).  With  the  addition  of  the 
boundary  B(18),  and  the  knowledge  of  similar  boundaries  in  this  frequency  range, 
the  following  conditions  become  the  dominant  constraints  on 

The  region  near  -10  dB  at  -70“  on  the  open  loop  grid  of  the  Nichols  chart  in 
Figure  4.4  is  the  focal  point.  To  meet  specifications  in  this  area;  the  magnitude  of 
■^p(jlO)  must  be  greater  than  the  magnitude  of  B(10),  but  there  is  also  the  bound 
at  u.’  =  18  rad/sec  requiring  lLc,(jT8)'  <  iB(18)|.  For  an  open  loop  phaise  angle  of 
approximately  -70'  the  rate  at  which  £y(ju.')  decreases  in  the  frequency  range  is 
rather  small.  This  results  in  a  nearly  horizontal  segment  of  Lo(j(^)  on  the  Nichols 
chart  after  suitable  additional  phase  angle  is  contributed  by  G(s).  If  additional 
phase  lag  is  added  at  frequencies  too  near  18  rad/sec  stability  bounds  similar  to 
J5(18),  unmarked  in  Figure  4.4,  are  violated.  For  this  reason  the  additional  phase 
lag  to  be  added  by  G(s),  and  required  to  decrease  the  magnitud*-  Lo(iw),  is 
delayed  until  higher  frequencies  for  a  given  complexity  of  G(5).  Th':  results  in  a 
higher  bandwidth  of  G{s). 

All  designed  G(s)  obtained,  using  the  above  reasoning,  are  felt  to  have  ex¬ 
cessive  bandwidth  and  gain.  This  problem  suggests  that  bounds  on  G{ju>),  rather 
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than  Lg{ju),  give  more  insight  into  the  problem  at  hand.  Bounds  on  G{ju)  are 
obtained  by  dividing  out  the  nominal  plant  Po{ju})  from  the  bounds  on  the  nomi¬ 
nal  loop  transmission  Lo{joj),  that  is  G(ju>)  =  ,  and  G(s)  is  designed  rather 

than  Xo(s).  However,  the  decision  is  made  to  continue  shaping  Lo(s}  in  the  fol¬ 
lowing  manner.  The  major  causes  of  excessive  bandwidth  and  gain  of  (^(s)  are  the 
modified  boundaries  B(8)  and  B(10)  of  Figure  4.4.  The  normal  tracking  bounds 
at  these  frequencies  for  ZLo(ji^)  >  —100°  are  well  below  -40  dB.  The  decision  is 
therefore  made  to  intentionally  violate  these  two  modified  bounds  (obtained  by  the 
method  in  the  Modified  Tracking  Bounds  subsection  of  this  chapter).  The  tracking 
boundary  at  B(0.5)  now  becomes  dominant  along  with  the  previously  mentioned 
stability  bounds,  allowing  for  a  much  smaller  gain  and  bandwidth  required  of  <j{5). 


The  above  reasoning  leads  to 


-2.6303  X  10^(5  +  5)2(5  +  40) 
s{s  +  50)  (s2  +  705  +  4900) 


(4.6) 


with  the  frequency  response  of  Figure  4.7,  which  results  in  a  very  modest  bandwidth 
of  Z/<,(s)  of  Wc  «  1.5  as  shown  in  Figure  4.8.  A  more  economical  compensator  can 
be  chosen  at  the  cost  of  greater  complexity. 


The  Design  of  the  Prefilter  F[s).  The  bounds  on  the  prenlter  F{s)  of  Fig¬ 
ure  2.1  can  be  generated  from  the  plant  templates  and  Nichols  chart  [13j;  or  by 
simulating  each  of  the  15  plants,  as  and  storing  the  maximum  and  mini¬ 

mum  magnitudes  at  each  frequency  as  the  variation.  The  second  technique  is  used 
in  obtaining  Figure  4.9.  The  constraints  on  F(5)  are  the  range  of  magnitudes,  at 
each  frequency,  which  ensures  that  the  response  of  is  within  the  bounds 

of  Figure  4.1.  The  resulting  bounds  on  F{s),  as  we.  as  the  first  F{s)  chosen. 


25(s  +  4) 

52  _  105  -  100 


(4.7) 


are  shown  in  Figure  4.10. 
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It  is  obvious  that  satisfying  the  bounds  on  F(s)  for  all  frequencies  results 
in  an  impractical  prefilter;  due  to  the  fact  that  has  a  greater  number  of 

excess  poles  than  the  two  of  the  lower  frequency  bound  of  Eq  (4.3).  The  required 
amplification  around  10  rad/sec  may  also  cause  problems  for  certain  cases.  One 
solution  is  to  obtain  another  lower  bound  in  the  frequency  domain  with  the  same, 
or  greater,  number  of  excess  poles  than  and  the  proper  time  response,  but 

this  extra  effort  is  not  necessary.  It  is  noted  that  this  could  have  been  done  at  the 
onset  of  the  problem.  A  simple  tradeoff  may  be  made  between  boundary  violation 
on  F(s)  and  obtaining  a  practical  prefilter.  This  tradeoff  does  not  present  a  real 
problem  since  the  violation  of  the  boundaries  can  be  deferred  in  frequency  until 
the  magnitude  of  L{jui)  is  very  small. 

The  MATRDCx  linear  simulation  of  the  final  design  is  presented  in  Chapter  V 
along  with  the  reasons  behind  modifying  the  prefilter.  Acceptable  time  response 
is  the  final  objective,  therefore  one  technique  is  to  use  the  slowest  prefilter  that 
results  in  acceptable  responses.  The  following  prefilter  is  chosen; 


F{s) 


SK, 

s  +  8 


(4.8) 


where  if/  is  a  nonlinear  gain  function  experimentally  determined  from  nonlinear 
simulations  (see  .4ppendi.x  C).  This  variable  gain  is  used  to  minimize  the  step 
error  coefficients  of  the  nonlinear  simulations  over  a  large  range  of  C*  commands, 
allowing  for  a  more  “linear"  response  for  the  system.  The  value  of  if/  =  1  is  used 
in  all  the  linear  simulations  of  the  inner  loop  SAS  on  MATRDCx. 


Summary 

This  chapter  presents  generation  of  the  bounds  on  Lo(s),  and  its  shaping  on 
the  Nichols  chart.  The  modification  to  the  usual  tracking  bounds  is  described  in 
detail  and  a  graphical  example  presented.  Two  choices  of  G{s)  and  their  design  are 
outlined,  as  well  as  the  logic  enabling  the  shortcomings  of  the  first  design  to  give 
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needed  insight  for  the  accepted  second  design.  Finally,  the  prefilter  is  designed  and 
a  discussion  on  practical  prefilters  is  presented. 
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V.  Simulation  of  the  Inner  Loop  SAS 


I 


Introduction 

This  chapter  presents  the  linear  and  nonlinear  simulations  of  the  stability 
augmentation  system  (SAS)  designed  in  Chapter  IV.  The  design  only  guarantees 
satisfaction  of  responses  for  Ig  and  2g  C*  commands  at  the  nominal  flight  con¬ 
dition  of  0.9M  at  20,000  ft;  however,  the  inherent  robustness  of  the  technique  is 
demonstrated  by  signiflcantly  varying  the  initial  conditions.  A  demonstration  of 
the  linearizing  effect  of  the  feedback  is  included.  The  inner  loop  is  then  modeled  as 
an  equivalent  LTI  plant  for  use  in  the  additional  pilot  compensation  of  Chapter  VI. 


MATRIXx  Linear  Simulation 


The  response  of  all  15  LTI  equivalent  plants  are  shown  in  Figure  5.1  using 
the  compensator  G(s)  of  Eq  (4.6)  and  the  prefilter 


F{s) 


25(5  +  4) 

5*  -r  10s  +  100 


of  Figure  4.8. 


The  MATRIXx  linear  simulations  are  used  only  as  a  design  aid  while  a  sat¬ 
isfactory  nonlinear  response  is  the  final  goal.  Therefore  further  inspection  of  the 
individual  responses  is  required.  In  the  proposed  design  of  G(s)  outlined  in  Chap¬ 
ter  IV,  in  which  a  limiting  cycle  of  the  nonlinear  YF-16  is  excited,  it  is  found 
that  the  equivalent  plant  #15  best  describes  the  nonlinear  aircraft  at  this  oper¬ 
ating  point.  This  conclusion  is  made  based  on  the  similarities  in  the  responses  of 
plant  #15  and  the  nonlinear  simulation  for  this  first  design  of  G(s).  The  response 
of  plant  #15  has  a  lightly  damped  sinusoidal  component  while  the  nonlinear  sim¬ 
ulation’s  response  shows  no  reduction  in  the  oscillation’s  magnitude  within  50  sec. 
For  the  final  design  the  response  of  plant  #15  in  Figure  5.1  has  a  30  %  overshoot 


I 
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Figure  5.1.  MATRDCx  Simulation  of  SAS  Over  Range  of  Uncertainty 

and  a  slight  amount  of  oscillation  before  reaching  its  final  value,  as  seen  in  trace 
of  Figure  5.2.  From  previous  experience  with  the  nonlinear  simulation,  it  is  ex¬ 
pected  that  the  response  of  the  nonlinear  aircraft  would  be  worse.  The  decision  is 
made  to  slow  down  the  prefilter,  specifically  to  remove  the  amplification  around 
10  rad/sec  from  Figure  4.8.  The  implementation  of  the  final  prefilter,  gives 
the  response  of  trace  in  Figure  5,2  for  plant  ^^15.  The  slower  prefilter  causes 
several  of  the  responses  of  the  set  P  to  violate  the  bounds  on  the  output;  however 
the  nonlinear  results  axe  the  end  objective.  The  simulations  are  then  continued 
using  the  first  order  prefilter. 

Nonlinear  Simulation 

The  compensation  is  then  implemented  in  Fortran  in  the  form  of  Figure  2.1 
page  2-2  where  the  output  of  G(s)  is  the  actuating  signal  to  the  elevator’s  actuators 
and  the  feedback  signal  is  C*.  The  initial  flight  condition  is  specified  and  the 
program  calculates  the  required  thrust  as  well  as  the  deflections  of  the  elevators 
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Figure  5.2.  Response  of  Plant  15  for  the  two  Prefilters 


o 


and  leading  edge  flaps  to  maintain  this  condition.  Once  the  equilibrium  conditions 
are  established,  the  internal  states  of  the  compensators  are  solved  for  their  trim 
values.  This  is  automatically  completed  prior  to  t  =  0.  The  simulation  is  then  run 
where  commands  are  given  as  perturbations  from  the  trim  values  and  outputs  are 
in  full  state  values. 

Nominal  Flight  Condition  Simulation.  The  first  run  is  a  Ig  commanded  C*  at 
the  nominal  flight  point  of  0.9M  and  20,000  ft,  as  shown  in  Figure  5.3.  Comparison 
of  this  nonlinear  response  with  the  MATRDCx  linear  response  of  plant  #15  for  the 
same  prefilter  validates  the  use  of  the  first  order  prefilter  for  this  case.  It  is  yet  to 
be  determined  if  the  slower  prefilter  adversely  affects  other  operating  points. 

The  design  guarantees  acceptable  responses  for  Ig  to  2g  commands  at  the 
nominal  flight  condition  and  Figure  5.4  shows  this  to  be  the  case.  The  responses 
are  normalized  to  the  commanded  input  over  the  range  of  0.5g  to  2.5g  steps  at 
increments  of  O.lg. 


The  nonlinear  gain  factor,  Kf,  mentioned  in  Chapter  IV  enables  the  responses 
of  Figure  5.4  to  be  normalized  to  commanded  input  rather  than  final  value.  Prior 
to  the  addition  of  this  gain,  the  steady  state  error  varied  nonlinearly  from  one 
command  to  another. 

It  is  desired  to  represent  the  inner  loop  SAS  as  a  single  linear  plant  at  the 
nominal  flight  condition  for  synthesis  of  the  outer  loop  of  Chapter  VI.  Although 
the  addition  of  the  feedback  greatly  linearizes  the  aircraft  response,  the  additional 
gain  increases  the  validity  of  this  approximation.  The  gain  Kf  is  scheduled  as  a 
piece-wise  linear  function  obtained  by  fitting  data  points  with  straight  line  approx¬ 
imations.  Each  data  point  is  equal  to  the  magnitude  of  the  command  divided  by 
the  steady  state  value.  The  function  is  specifically  derived  for  0.9M  at  20,000  ft 
and  is  contained  in  Appendix  C. 

Expanding  the  Envelope.  At  the  onset  of  the  design,  equivalent  LTI  plajits 
are  only  obtained  at  the  one  nominal  flight  condition  for  Ig  and  2g  responses. 
Therefore,  the  design  is  only  guaranteed  to  work  over  this  range.  To  guarantee  a 
larger  flight  envelope  requires  obtaining  enough  equivalent  plants  at  various  flight 
conditions  to  fully  represent  the  nonlinear  aircraft  over  that  region  of  operation. 
However,  the  current  set  of  plants  P  models  much  of  the  nonlinear  characteristics 
already  and  may  in  fact  be  sufficient.  This  statement  is  further  supported  by 
the  results  of  several  hundred  simulations  consisting  of  step  responses  of  various 
magnitudes  and  flight  conditions. 

The  simulation  envelope  spans  the  region  of  5,000  to  30,000  ft  with  speeds 
from  0.5  to  0.9M  and  commands  up  to  9g’s,  all  simulations  shown  start  trimmed 
for  straight  and  level  flight.  Figure  5.5  shows  the  results  for  commands  of  Ig  to  9g’s 
normalized  to  final  value,  while  representative  results  at  25  flight  conditions  are 
shown  in  Appendix  C.  The  response  for  the  mtiximum  command  at  the  nominal 
point  is  shown  in  Figures  5.6  and  5.7  where  all  angles  are  in  degrees,  rates  are  in 
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Figure  5.5.  Normalized  C*  Responses  for  l-9g  Commands  at  Nominal 
deg/sec,  altitude  in  feet,  and  true  air  speed  in  ft /sec. 

The  decision  as  to  whether  or  not  a  particular  response  at  each  flight  con¬ 
dition  is  included  in  the  envelope  is  made  by  checking  the  normalized  response 
against  the  specified  tracking  bounds  of  Figure  3.2  page  3-4.  This  envelope  is 
shown  in  Figure  5.8  with  a  sampling  of  these  limits  in  Table  5.1.  Each  loci  of 
points  represents  the  maximum  command  resulting  in  an  acceptable  response  for 
a  given  altitude  as  airspeed  is  varied.  It  is  felt  that  the  major  constraint  on  the 
operating  region  is  due  to  aerodynamics  and  not  the  controller.  The  decrezising 
maneuverability  of  the  aircraft  with  dynamic  pressure,  along  with  the  decrease  in 
available  thrust  at  higher  altitudes,  explain  the  results  well. 

Two  comments  on  generating  the  acceptable  response  envelope  require  clar¬ 
ification.  First,  the  responses  are  normalized  to  their  final  value  not  commanded 
value.  Therefore,  the  envelope  does  not  suggest  that  the  commanded  C*  can  be 
maintained  at  every  flight  condition.  Second,  several  responses  at  0.5M  violated 
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Figure  5.6.  Surface  Deflections  and  Rates  for  9g  Command  at  Nominal 


Altitude  (ft) 

Mach  Number 

0.5  0.6  0.7  0.8  0.9 

5,000 

7.3 

9.0 

9.0 

9.0 

9.0 

10,000 

5.9 

9.0 

9.0 

9.0 

9.0 

15,000 

■ai 

6.5 

9.0 

9.0 

9.0 

20,000 

3.7 

5.0 

9.0 

9.0 

30,000 

2.0 

3.3 

4.2 

4.2 

— 1 

Table  5.1.  Maximum  Commands  in  g’s  for  Acceptable  Responses 
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ACCEPTABLE  AESPONSE  ENUELOPC  5*000  TO  30.000  FEET 

Figure  5.8.  Flight  Envelope  of  Acceptable  SAS  Responses 

the  bounds  for  small  commands,  then  achieved  acceptable  responses  for  larger 
command  up  to  the  point  of  the  maximum  command  listed  in  the  envelope  (see 
Appendix  C) .  This  is  due  to  the  slow  prefilter  chosen  and  can  be  eliminated  by  in¬ 
creasing  the  bandwidth  of  the  prefilter.  It  should  be  noted  that  the  slow  responses 
occur  while  expanding  the  envelope  and  axe  not  characteristic  of  the  nominal  flight 
condition. 

Equivalent  Inner  Loop 

As  mentioned  in  the  discussion  of  the  nonlinear  prefilter  gain,  it  is  desired  to 
represent  the  designed  SAS  as  a  new  equivalent  linear  plant.  This  simplifies  the 
synthesis  of  the  additional  pilot  compensation  of  Chapter  VI. 

The  program  used  to  generate  the  set  P  is  used  to  obtain  the  equivalent  plant 


P.(s)  = 


97.5S70(s  -f  0.8703)(s2  +  3.9438  x  yiO.3674) 

(s  ^  0.8352)  (s-  +  6.3252  ±  y8.8315)(s2  +  8.7786  ±  i5.4357) 
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A* 


Figure  5.9.  Outputs  of  Pe{s)  and  YF-16  for  Step  Input 


The  time  histories  used  are;  the  C*  response  for  a  commanded  Ig  at  the  nominal 
flight  condition,  while  the  input  is  a  unit  step  instead  of  the  signals  to  the  actuators. 
The  response  of  both  the  nonlinear  simulator  and  the  equivalent  SAS  plant  to  a 
Ig  step  axe  shown  along  with  the  error  in  Figure  5.9. 


Summary 

This  chapter  presents  the  MATRDCx  linear  simulation  of  the  SAS  designed  in 
Chapter  IV  as  well  as  the  reasons  for  modifying  it.  The  modification  is  validated  by 
comparing  the  resulting  linear  simulation  with  the  nonlinear  response.  A  nonlinear 
prefilter  is  then  presented  to  better  linearize  the  inner  loop  final  values  for  use  in 
the  synthesis  of  Chapter  VI.  The  robustness  of  the  inner  loop  is  demonstrated 
by  nonlinear  simulations  that  drastically  expand  the  flight  envelope.  Finally,  an 
equivalent  LTI  SAS  plant  is  chosen  to  represent  the  design  up  to  this  point. 
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VI.  Pilot  Compensation 


Introduction 

This  chapter  presents  the  pilot  compensation  to  be  used  in  conjunction  with 
the  SAS  designed  in  Chapter  IV.  The  synthesis  of  this  compensation  is  based 
upon  obtaining  a  desirable  response  from  the  aircraft  while  reducing  the  workload 
of  the  pilot.  The  technique  used  is  based  upon  the  Neal-Smith  pilot  model  for 
compensatory  tracking.  The  compensation  is  designed  on  MATRIXx  followed  by 
nonlinear  simulations  in  Fortran.  A  fourth  order  Fade  approximation  of  a  pure 
time  delay  is  used  in  the  Fortran  simulations  and  the  effects  of  this  approximation 
are  shown. 

The  Pilot  Model 

([.•  The  form  of  the  mathematical  model  used  in  the  design  of  the  remaining 

compensation  is  shown  in  Figure  6.1.  The  equivalent  plant  Pe(5)  is  that  of  Eq  (5.1) 
representing  the  SAS  designed  in  Chapter  IV.  The  model  used  to  represent  the  pilot 
is  that  used  in  the  technical  report  by  Peter  Neal  and  Rodger  Smith  [18]  to  predict 
pilot  ratings  for  a  longitudinal  compensatory  tracking  task.  The  configuration  of 
Figure  6.1  allows  only  a  single  degree  of  freedom,  Fp{s),  to  achieve  the  desired 
response. 

The  dynamic  model  representing  the  pilot  is  taken  directly  from  the  Neal- 
Smith  report  [18|  as: 

The  pilot  is  represented  by  a  variable  gain  Kp,  a  time  delay,  and  a  variable  first 
order  compensation  network.  The  factors  Vp^  and  Tp^  represent  the  adaptive  com¬ 
pensation  provided  by  the  pilot  to  accomplish  a  specific  task.  The  time  delay 
includes  the  time  required  for  the  pilot  to  sense  a  change  in  C*,  decide  what  action 


(6.1) 
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Figure  6.1.  Form  of  Pilot  Compensation 

to  take,  and  the  neuromuscular  lags  (I8j.  The  time  delay  of  0.3  seconds  used  may 
be  excessive  since  the  YF-16  uses  a  force  stick  and  the  displacement  of  the  pilot’s 
arm  is  negligible.  The  neuromuscular  lags  could  therefore  be  neglected,  but  the 
objective  is  to  present  a  synthesis  technique  so  the  value  of  0.3  seconds  is  used 
throughout  this  thesis.  The  purely  mental  time  delay  of  a  human  is  closer  to  0.2 
seconds  ilS.. 

The  block  diagram  of  Figure  6.1  is  considered  a  compensatory  model  in 
that  the  pilot  operates  only  on  the  difference  of  commanded  and  actual  C*.  It  is 
obvious  that  in  real  life  this  is  not  the  case,  the  pilot  also  uses  many  other  inputs 
in  determining  the  control  that  he  applies.  However,  this  model  has  been  shown  to 
adequately  explain  the  more  important  aspects  of  tracking  [18],  and  is  relatively 
simple  to  analyze. 
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Transforming  Analysis  to  Synthesis 

The  work  in  pilot  modeling  to  date  has  been  used  to  predict  the  pilot  rating 
of  a  designed  flight  controller;  if  unacceptable  results  were  predicted,  modiflcations 
could  be  made  to  the  prefilter  (l|.  The  objective  of  this  thesis  is  to  take  into  ac¬ 
count  the  prediction  of  pilot  rating,  given  by  the  Neal-Smith  criterion,  and  provide 
compensation  a  priori  to  give  an  acceptable  predicted  response. 

The  goal  of  this  compensation  is  to  minimize  the  additional  compensation,  or 
workload,  required  of  the  pilot  to  meet  good  closed  loop  performance  requirements. 
The  time  delay  in  Eq  (6.1)  is  inherent  to  the  pilot,  it  is  the  gain  Kp  and  the  first 
order  compensation  network  that  must  be  generated  by  the  pilot.  Minimizing  this 
compensation  would  enable  the  pilot  to  act  only  as  a  time  delay  giving,  Kp  =  \  and 
Tp,  =  Tpj  =  0.  The  proposed  synthesis  technique  is  to  perform  a  prediction  zmalysis 
similar  to  that  of  Neal-Smith,  while  forcing  an  acceptable  rating.  The  resulting 
pilot  model  may  be  of  greater  complexity  than  Eq  (6.1)  but  must  contain  the  time 
delay  term.  The  resulting  model  is  then  equal  to  Pilot  •  Fp[s)  of  Figure  6,1.  Since 
it  is  desired  that  the  pilot  provide  no  additional  lead-lag  compensation,  Fp[s)  is 
built  into  the  control  system  and  is  set  equal  to  the  resulting  pilot  model  obtained 
above  excluding  the  time  delay. 

Synthesis  of  the  Pilot  Compensation  Fp[s) 

As  outlined  earlier,  the  design  of  Fp{s)  is  based  upon  conducting  a  prediction 
of  pilot  rating  using  a  modified  Neal-Smith  criterion.  The  modification  is  that 
the  designer  is  free  to  model  the  pilot  with  any  lead  and/or  lag  compensation, 
instead  of  the  ‘'optimuni  ccmpexisalion  presented  in  the  technical  report  [18]. 
This  freedom  allows  one  to  obtain  an  acceptable  prediction  in  cases  which  would 
be  impossible  with  Eq  (6.1).  This  freedom  must  be  weighted  against  achieving  a 
practical  compensator  Ep(s). 
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Pilot’s  View  of  Good  Tracking  Performance.  The  Neal-Smith  report  [18] 
gives  an  excellent  explanation  of  converting  “what  the  pilot  sees”  into  usable  track¬ 
ing  performance  standards,  to  which  the  reader  is  referred.  Sections  6.2  and  6.3  of 
the  report  are  summarized  below  by  presenting  the  performance  standards  along 
with  definitions  used  in  the  technique  [18].  All  references  to  the  closed  loop  are  for 
of  Figure  6.1. 

Bandwidth  (wbh?);  Bandwidth  is  defined  as  the  frequency  for  which  the  phase 
angle  of  the  closed  loop  is  equal  to  —90°.  It  is  a  measure  how  quickly  the 
pilot  can  move  the  aircraft’s  nose  towards  the  target. 

Droop:  Droop  is  defined  as  the  maximum  excursion  of  the  closed  loop  magni¬ 
tude  below  the  0  dB  line  for  frequencies  less  than  ubw-  In  absence  of  large 
oscillations,  droop  is  a  measure  of  how  slowly  the  aircraft  settles  down  on 
target. 

Standards  of  Performance:  A  minimum  bandwidth  ubw  of  3.5  rad/sec  for  air¬ 
speeds  above  those  causing  wing  buffet  and  a  maximum  droop  of  3  dB. 

PIO  Tendency:  The  tendency  for  pilot  induced  oscillations  is  defined  in  terms  of 
the  Bode  magnitude  of  any  closed  loop  resonant  peak  that  results  from  the 
pilot's  efforts  to  achieve  the  performance  standards. 

Pilot  Compensation:  The  pilot’s  physical  and  mental  workload  required  to  achieve 
the  standard  is  defined  in  terms  of  the  phase  angle  of  his  additional  compen¬ 
sation  at  w  =  wgiv. 

A  visualization  of  the  closed  loop  standards  on  bandwidth,  droop  along  with 
a  PIO  limit  of  3  dB  are  shown  as  closed  loop  contours  on  the  Nichols  chart  of 
Figure  6.2. 

In  addition  to  the  above  standards  of  performance,  this  thesis  proposes  fur¬ 
ther  constraints  on  the  response.  These  modifications  to  the  standards  are  gener- 
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Open  Loop  Magnitude  (dB) 


c« 


ated  from  analyzing  the  pilot  comments  of  the  reference  report  [18].  Their  com¬ 
ments  suggest  that  a  response  with  an  initial  overshoot  but  no  resulting  undershoot 
of  the  final  value  is  preferred.  Such  responses  were  consistently  preferable  to  first 
order  responses  or  those  that  oscillated  about  the  final  value  like  a  pure  second 
order  response.  The  amount  of  undershoot  is  related  to  the  amount  of  droop  in 
the  closed  loop  frequency  response.  With  the  above  reasoning  and  the  fact  that 
the  original  3  dB  droop  limit  is  chosen  “somewhat  arbitrarily”  [18]  it  is  decided  to 
minimize  the  droop  to  nearly  0  dB. 

Shaping  the  Response.  Since  it  is  desired  to  have  the  pilot  act  only  as  a 
time  delay,  the  pilot  model  of  Figure  6.1  is  merely  for  the  remainder  of  this 
analysis.  The  open  loop  transfer  function  for  the  technique  is  then  defined  as: 

Open  Loop  =  0(s)  =  (6.2) 

With  only  a  single  degree  of  freedom  the  poles,  zeros,  and  gain  of  Fp(5)  are  used  to 
shape  the  closed  loop  directly.  Using  the  new  specification  of  droop,  the  following 
is  an  outline  of  the  shaping  of  on  the  Nichols  chart  (see  Figure  6.2). 

The  response  should  “hug”  the  0  dB  M  contour  until  crossing  the  closed  loop 
phase  angle  contour  of  —90°  at  a  frequency  of  u;  =  ujbw  =  3.5  rad/sec.  Once  the 
standards  are  meet  for  w  <  3.5  rad/sec  the  PIO  tendency  is  the  dominant  factor. 
The  Neal-Smith  Criterion  gives  the  maximum  closed  loop  resonance  resulting  in 
a  good  rating  as  3  dB.  Therefore  the  3  dB  M  contour  is  the  dominant  factor  for 
u!  >  3.5  rad/sec. 

The  pilot  model  and  Pe(s)  are  fixed,  thus  a  starting  point  in  the  design  is 
the  plot  of  the  open  loop  response  on  the  Nichols  chart  for  Fp{s)  —  1.  However,the 
desired  closed  loop  is  to  hug  the  0  dB  contour;  therefore,  the  first  try  is  to  place  a 
pole  at  the  origin.  The  response  for  Fp{s)  =  Ms  shown  in  Figure  6.3.  This  puts 
the  open  loop  0{ju)  for  u;  =  3.5  rad/sec  at  -12.1  dB  and  —192°.  The  compensator 
must  contribute  approximately  55°  at  this  frequency. 
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Figure  6.3.  Starting  Point  for  the  Design  of  Fp(s) 


The  problems  presented  by  the  time  delay  are  obvious  in  the  figure.  The 
characteristics  of  being  0.3a/  rad.  Having  a  gain  of  0  dB  at  all  frequencies 
but  a  phase  angle  contribution  of  —0.3  x  =  —17.19  deg/rad,  for  example 
—1.719“  at  0/  =  0.1  rad/sec  and  —171.9“  at  a/  =  10  rad/sec.  These  characteristics 
cause  the  open  loop  magnitude  \0{ju))\  to  be  nearly  horizontal  for  all  frequencies 
above  approximately  20  rad/sec.  However,  there  is  still  room  for  improvement  over 
the  pilot’s  baindwidth  of  0.5  to  10  rad/sec  [18].  This  near  0  dB/dec  slope  of  the 
open  loop  can  cause  extreme  gain  sensitivity  and  high  frequency  resonant  modes. 

One  example  depicting  both  of  these  conditions  is  shown  in  Figure  6.4  which 
is  taken  directly  from  the  Neal-Smith  report.  This  figure  shows  the  open  loop  0{juj) 
frequency  response  for  0.6  <  a/  <«  21  rad/sec  on  the  Nichols  chart.  The  solid  line 
is  for  0.6  <  0/  <  10  rad/sec  and  the  dotted  line  is  for  15  <  a/  <«  21  rad/sec.  The 
solid  line  shows  the  response  over  the  pilot’s  bandwidth  with  a  resonance  of  1  dB 
at  approximately  6.5  rad/sec.  An  increase  of  only  4  dB  in  the  open  loop  produces 


Open  Loop  Magnitude  (dB) 


a  resonance  of  12  dB,  representing  a  very  lightly  damped  response.  Visible  in 
this  frequency  range  is  that  the  slope  of  the  open  loop  function  0{ju>)  is  positive 
between  7  and  15  rad/sec  causing  the  high  frequency  resonance.  This  problem 
is  seen  by  observing  several  points  on  the  response  of  0{ju)  in  Figure  6.4.  The 
curve  shows  the  closed  loop  magnitudes  for  w,  «  10,15,17,21  rad/sec  to  be  ap¬ 
proximately  —5,  —6,  -2,-13  dB  respectively.  The  resonant  peak  in  this  frequency 
range  represents  a  poorly  damped  mode  at  u;  «  17  rad/sec  and  can  be  seen  in  a 
Bode  plot  of  the  closed  loop  magnitudes  for  the  above  five  frequencies.  The  pilot 
comments  on  the  handling  of  the  flight  control  system  represented  in  Figure  6.4  of 
“difficulties  in  acquiring  a  target”  and  “oscillations  about  the  target”  are  in  agree¬ 
ment  with  the  analysis.  Therefore  it  is  concluded  that  the  slope  of  the  open  loop 
transfer  function  is  to  be  kept  negative  until  the  magnitude  is  sufficiently  small  to 
avoid  the  above  problems.  However  the  phase  required  of  the  compensator  may  be 
so  large  that  this  is  not  possible,  for  such  cases  the  slope  is  kept  negative  as  long 
as  possible  while  meeting  low  frequency  specs. 


In  the  attempt  to  minimize  droop  and  PIO  tendencies,  while  achieving  some 
overshoot  and  a  bandwidth  of  3.5  rad/sec,  the  following  compensator  is  chosen: 

761.5615(5^  +  16.65  +  36) 


F  fsl  = 

’  5(52  -  1005  -  1002) 

The  system  frequency  response  is  shown  in  Figures  6.5  and  6.6. 


(6.3) 


Simulation  of  the  Pilot  Compensation 

The  MATRIXx  linear  time  responses  are  used  simultaneously  with  the  Nichols 
chart  in  obtaining  the  final  design.  The  step  response  is  as  desired,  an  initial  20  % 
overshoot  with  minimal  undershoot.  To  acquire  zero  undershoot  with  a  20  %  over¬ 
shoot  requires  a  reduction  in  bandwidth  which  is  considered  more  undesirable  than 
the  undershoot.  This  is  due  to  the  amount  of  lead  phase  angle  required  to  achieve 
a  bandwidth  uibw  —  3.5  rad/set  while  maintaining  a  reasonable  rate  of  decrease 
in  the  magnitude  of  0{jijj).  The  MATRIXx  linear  simulations  using  a  pure  delay 
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Figure  6.5.  Response  and  Standards  of  Performance  for  Final  Design 


Figure  6.6.  Closed  Loop  Frequency  Response 
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Figure  6.7.  Effect  of  Using  Fade  Approximation 


for  the  pilot  model  show  zero  output  for  the  first  0.3  seconds  with  time  specs  of 
tp  =  0.8  sec  and  =  1.6  sec. 


Again  the  true  test  is  that  of  the  nonlinear  simulation.  Instead  of  writing 
the  Fortran  code  to  approximate  the  time  delay  it  is  decided  to  use  a  Fade  ap¬ 
proximation  of  the  function  e'*  (see  Appendix  D).  Figure  6.7  shows  the  MATRIXx 
simulation  of  the  design  using  both  the  time  delay  and  a  fourth  order  Fade'  ap¬ 
proximation  along  with  the  error  caused  by  using  the  approximation.  Most  of  the 
error  is  during  the  first  0.3  seconds  where  the  approximation  has  a  nonminimum 
phase  response  while  the  delay  hais  zero  output.  The  approximation  is  judged  to 
be  sufficient  and  the  nonlinear  simulations  begun. 


The  first  nonlineax  simulation  is  for  a  Ig  step  commanded  initiated  from 
straight  and  level  flight  at  0.9M  and  20,000  ft.  The  response  is  shown  in  Figure  6.8 
and  is  judged  to  be  excellent,  again  20  %  overshoot  with  negligible  undershoot  and 
a  small  rise  time  after  the  initial  time  delay.  The  effects  of  the  approximation’s 
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Figure  6.8.  Nonlinear  Response  to  =  Ig 

nonminimum  phase  response  on  the  aerodynamics  are  considered  negligible  due 
to  their  magnitude  and  duration.  It  is  again  noted  that  this  nonminimum  phase 
characteristic  should  be  absent  when  a  genuine  time  delay  is  used. 

Next  step  is  to  answer  the  question  of  how  to  accomplish  other  simulations 
without  voiding  the  standards  of  performance  or  pilot  model  proposed  by  the  Neal- 
Smith  report.  Any  pilot  modeling  must  be  constrained  to  a  specific  task  [18].  The 
data  for  the  Neal-Smith  analysis  was  taken  during  fine  tracking  tasks,  not  for  large 
acquisition  commands.  The  decision  is  made  to  simulate  Ig  steps  while  varying 
the  initial  trim  conditions.  Additional  simulations  starting  from  coordinated  turns 
trimmed  to  load  factors  of  2  to  5g’s  are  found  to  give  similar  results  as  shown  in 
Figure  6.9. 

The  question  remains  as  to  the  sensitivity  of  the  system  to  changes  in  the 
bandwidth  while  performing  the  same  task.  That  is,  what  are  the  effects  of  this 
additional  compensation  when  the  pilot  flies  to  achieve  a  bandwidth  other  than 


6-12 


time 


Figure  6.9.  Simulations  for  Trimmed  Load  Factors  of  Ig  to  5gs 

3.5  rad/sec.  This  is  accomplished  by  referring  to  Figure  6.5.  If  a  reduction  in 
bandwidth  is  desired,  the  pilot  need  only  lower  his  gain,  this  has  a  stabilizing 
effect. 

To  increase  the  bandwidth  the  pilot  must  use  gain,  lead  compensation,  or 
a  combination  of  the  two.  Pilots  seem  to  have  no  problem  with  generating  up 
to  25°  of  lead  compensation  and  feel  comfortable  with  doing  so  [18].  Without 
the  compensation  ip(s)  the  pilot  may  have  used  all  this  lead  compensation  and 
more,  although  the  handling  may  have  been  judged  good  using  the  Neal-Smith 
Criterion.  However  this  is  based  upon  achieving  a  bandwidth  of  3.5  rad/sec.  To 
increase  this  bandwidth  the  pilot  would  have  to  use  mostly  gain,  if  he  had  used 
up  his  lead  generating  ability,  which  gives  the  gain  sensitivity  problem  disc\issed 
with  Figure  6.4.  The  addition  of  Fp{s)  allows  the  pilot  to  achieve  a  bandwidth 
of  3.5  rad/sec  without  generating  any  additional  lead.  This  comfortable  region  of 
additional  lead  is  then  used  to  increase  the  bandwidth  of  the  entire  closed  loop 
system  avoiding  the  gain  sensitivity  problem. 


Summary 


This  chapter  presents  the  design  of  the  outer  loop  of  Figure  6.1,  consisting 
of  pilot  compensation  based  upon  the  Neal-Smith  Criteria.  The  objective  of  min¬ 
imizing  the  compensation  required  of  the  pilot  is  accomplished  by  providing  the 
compensation,  to  achieve  a  bandwidth  of  3.5  rad/sec  with  a  good  response,  with 
Fp{s).  This  allows  the  pilot  to  act  as  a  pure  time  delay.  The  nonlinear  simulations 
presented,  using  a  Fade  approximation  of  the  time  delay,  show  excellent  results  for 
initial  load  factors  of  up  to  5gs.  Finally,  a  discussion  of  the  effects  of  bandwidth 
show  that  the  closed  loop  bandwidth  for  the  tracking  task  has  been  increased. 
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VII .  Conclusions  and  Recommendations 


Discussion 

This  thesis  demonstrates  the  application  of  Nonlinear  Quantitative  Feedback 
Theory  to  flight  control  design  and  presents  a  compensation  design  technique  which 
reduces  pilot  work  load.  The  nonlinear  QFT  technique  is  applied  to  the  YF-16 
aircraft  which  results  in  an  inner  loop  SAS  to  control  the  output  C*  and  guarantees 
that  the  response  lies  within  the  bounds  suggested  by  Tobie  [20].  However,  this 
inner  loop  design  does  not  include  the  characteristics  of  the  human  pilot  a  priori, 
and  such  is  the  case  with  nearly  all  current  preliminary  flight  control  designs  [ij. 
The  characteristics  of  the  human  pilot  are  usually  not  dealt  with  until  the  man- 
in-the-loop  simulation  phase  of  final  control  system  design.  Waiting  until  this 
late  point  in  the  system’s  development  to  attack  this  problem  is  extremely  cost 
inefficient  in  two  ways.  First,  due  to  the  high  cost  of  such  simulations  delay  must  be 
keep  to  a  minimum.  Secondly,  modifications  to  a  system  so  far  into  it’s  development 
are  very  costly.  Ignoring  the  inherent  time  delay  of  the  human  pilot  results  in 
extensive  adjustments  to  the  original  design  to  achieve  satisfactory  pilot  ratings. 

Instead  of  fine  tuning  the  design  during  this  expensive  phase  of  system  design, 
an  adjustment  is  made  prior  to  these  man-in-the-loop  simulations  with  additional 
pilot  compensation.  This  pilot  compensation  allows  the  controller  to  achieve  good 
handling  qualities  predictions  for  simulations  including  the  inherent  time  delay  of 
the  human  pilot. 

Extensive  use  of  CAD  tools,  not  included  in  this  thesis,  greatly  reduces  the 
calculations  and  book  keeping  required  in  the  synthesis  process.  Listings  of  the 
Fortran  code  and  MATRIXx  routines  may  be  made  available  through  the  following 
points  of  contact: 
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Convert  Time  Histories  to  Transfer  Functions:  Professor  Isaac  M.  Horowitz, 
Department  of  Electrical  Engineering  and  Computer  Sciences,  University  of 
California  at  Davis. 

MATRIXx  Macros:  Professor  Constantine  H  Houpis,  Electrical  and  Computer 
Science  Engineering  Department,  Air  Force  Institute  of  Technology,  WPAFB 
Dayton  Ohio. 


Conclusions 

Based  upon  the  work  performed  on  this  thesis  and  the  results  obtained,  the 
following  conclusions  are  made. 


1.  A  nonlinear  plant  can  be  represented  by  a  set  P  of  linear  time- invariant  plants 
for  use  in  a  design  problem.  The  method  of  Chapter  III  is  found  to  describe 
the  nonlinear  YF-16  rather  well. 


2. 


Nonlinear  Quantitative  Feedback  Theory  can  be  successfully  applied  to  the 
synthesis  of  flight  controllers  resulting  in  fixed  compensators  that  yield  ac¬ 
ceptable  responses  over  a  large  range  of  nonlinear  parameter  variation. 


3.  Linear  simulations  are  an  effective  design  aid  but  final  analysis  of  the  control 
of  the  nonlinear  plant  must  be  performed  using  the  nonlinear  simulations. 
The  linear  simulations  are  invaluable  for  the  insight  that  they  provide  but  this 
insight  can  only  be  gained  by  comparing  them  to  the  nonlinear  simulations. 

4.  Use  of  the  MATRIXx  routines  in  the  synthesis  process  greatly  reduces  the 
workload  and  calculation  errors,  however  they  are  written  for  this  specific 
thesis.  There  is  a  definite  need  for  a  generic  computer  aided  design  (CAD) 
package  to  ease  the  design  burden  associated  with  QFT  and  to  speed  the 
design  process. 

5.  It  is  possible  to  design  compensation  for  the  task  of  reducing  pilot  workload 
and  therefore  increasing  the  system  bandwidth.  The  technique  developed 
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in  this  thesis  will  work  for  other  tasks.  It  is  necessary  to  emphasize  that 
the  closed  loop  performance  requirements  are  task  specific,  and  the  expected 
performance  is  expected  to  vary  depending  on  the  particular  tzisk. 


Recommendations 

After  completion  of  the  of  the  stability  augmentation  system  design  for  the 
YF-16  and  the  pilot  compensation,  followed  by  analyzing  the  simulation  results, 
one  task  remains  to  be  completed.  The  technique  must  be  validated  by  actual 
man-in-the-loop  simulations  such  as  on  AFIT’s  SIMSTAR  hybrid  computer.  The 
real-time,  full  flight  envelope  simulator  is  nearing  completion  [2,16].  An  original 
objective  of  this  thesis  effort  was  to  include  such  simulations;  however,  upgrades  to 
the  system  have  prevented  it.  Beyond  this  immediate  suggested  task,  three  distinct 
areas  for  future  study  are  proposed. 


Extend  the  nonlinear  SISO  plant  to  that  of  the  nonlinear  MIMO  plant  with 
an  emphaisis  on  survivability  and  on  acceptable  responses  while  experienc¬ 
ing  various  surface  failure  conditions.  Such  work  is  underway  based  upon 
linearized  plant  models;  however,  using  the  nonlinear  plant  should  make  the 
results  more  realistic. 


2.  Expand  the  area  of  pilot  compensation.  Given  any  existing  flight  controller, 
use  the  developed  technique  for  a  variety  of  pilot  tasks,  where  each  task  has 
different  performance  requirements.  The  controller  can  be  adapted  for  a  given 
task  and  a  new  compensator  designed.  Included  in  this  research  should  be  the 
ability  to  perform  real-time  man-in-the-loop  simulations  such  as  on  AFIT’s 
SIMSTAR  F-16  simulator.  This  type  of  simulation  can  validate  the  technique 
proposed  in  this  thesis  by  comparing  handling  qualities  ratings  for  the  same 
aircraft  before  and  after  the  implementation  of  the  pilot  compensation  for 
different  tasks. 
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3.  Develope  a  generic  interactive  CAD  package  devoted  to  QFT.  Many  useful 
routines  have  been  developed  by  researchers  for  their  work;  however,  their 
objective  is  to  complete  a  specific  design.  Therefore  these  routines  are  written 
for  the  specific  problem  at  hand  with  no  time  to  check  out  other  possibilities 
that  can  arise  in  other  situations.  A  full  time  effort  devoted  solely  to  such  a 
CAD  package  will  be  a  step  in  the  right  direction. 

Summary 

This  chapter  reviews  the  design  and  results  of  this  thesis.  Recommendations 
are  made  for  future  projects  in  nonlinear  QFT,  pilot  compensation,  and  develop¬ 
ment  of  a  QFT  CAD  package. 
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Appendix  A.  Equivalent  LTI  Plants 


The  following  axe  the  LTI  equivalents  of  the  nonlinear  plant  response,  out¬ 
lined  in  Chapter  III. 

C‘  Transfer  Functions  for  Ig  Responses 
Plant  #1:  fit  with  2,4 

C*(s)  _  -7.8130(s  +  2.2038)(5  +  11.2499) 

es,{s)  ~  (5  -h  1.440l)(s  +  2.535  ±  j3.02Z2){s  +  54.9636) 

Plant  ff2:  fit  with  1,3 

C*(s)  _  -12.150l{s  + 4.4432) 

Cfe(s)  ”  (5 +  2.1664  ±>3.1576) (s  +  23.6823) 

Plant  #3:  fit  with  2,4 

C*(s)  _  -5.3139(5 +  3.6184)(5  + 4.4139) 

e«.(5)  ~  (5  +  1.6627  ±>2.5911)(5  +  5.7862  ±>4.8427) 

Plant  #4;  fit  with  2.4 

C*{s)  _  -5.0612(s  + 2.4154  ±>2.1172) 

eec{s)  ~  (5  +  1.6203  ±  >2.6521)(5  +  4.8324  ±  >3.3295) 

Plant  ffS:  fit  with  3,5 

C*(5)  _  -5.305(5  -  0.3668)(5  +  2.7128)(5  +  37.5146) 

ese{s)  ~  (5  ±0.3301)(5  ^  1.816  ±>2.9212)(5  + 9.2333  ±>14.9602) 

Plant  #6:  fit  with  3,5 

C*(5)  _  -8.4862(5  +  1.4992)(5  +  3.5374  ±  >9.092) 

eet{s)  ~  (5  +  1.1349)(5  +  3.158  ±  >3.5824)(5  -h  4.2219)(5  +  70.153) 
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Plant  t^7:  fit  with  4,6 

C*(5)  _  -7.557(5  +  0.4904)(5  +  4.88271  ±  jl7.0975)(5  +  6.3816) 

eseis)  ~  {3  +  0.4382) (5  +  2.4838  ±  j2.973){s  +  6.8833  ±  il3.8218)(s  +  29.3342) 

Plant  ^8:  fit  with  3,5 

C*(s)  -5.6663(5 +  0.2243)(s  +  2.4005)(s  + 27.2064) 

C«,(5)  “  (5  +  0.1912) (5  +  1.6578  ±  j2.9.563)(5  +  8.197  ±  j  12.3188) 

C*  Transfer  Functions  for  Sg  Responses 
Plant  #9:  fit  with  2,4 

C*(5)  _  -7.0856(5  +  1.3909) (5  +  15.8178) 

e«.(5)  “  (5  +  0.8592)  (5  +  2.8374  ±  j2.3273)(5  +  35.7299) 

Plant  #10:  fit  with  2,4 

C*(5)  _  -7.2503(5  +  3.587  ±j3.1188) 

e«e(s)  ■"  (5  +  1.5046  ±  y2.3305)(5  ^  4.7352  ±  ;5.923l) 

Plant  #11:  fit  with  3,5 

C*(5)  _  -7.3788(5  +  0.9628)(5  +  3.7751  ±  >6.8517) 

e«*(s)  ~  (5  +  0.653)(5  4-  3.0254)(5  +  3.8841  ±  >3.815)(5  +  19.3303) 

Plant  #12:  fit  with  2,4 

C‘(5)  _  -8.9094(5  -  0.977  =  >2.4645) 

€s,{s)  ~  (5  -  0.8337  X  >2.5676)(5  -f  3.0528)(5  -  7.6762) 

Plant  #13:  fit  with  3,5 

C*  (5)  _  -6.0796(5  +  1.1201  ±  >3.2202)  (5  +  9.4709) 

6^,(5)  ~  (5  +  0.9523)(5  +  1.0421  ±  >3.0456)(5  +  10.2039  ±  >8.8792) 

Plant  #14:  fit  with  3,5 

C*(5)  _  -8.8134(5  +  0.6449)(5  +  7.0771  ±  >6.0902) 

e«e(5)  ~  (5 +  0.4816)(5  + 2.8856  ^>2.2334)(5  + 9.6376  ±>10.1537) 

Plant  #15:  fit  with  4,6 

C*(5)  _  -10.164(5  +  0.3891)(5  +  1.946  ±>20.3273)(5  ±3.0283) 

e«e(5)  ~  (5  ±  0.3002)(5  ±  1.7901  ±  >2.1313)(5  ^  3.4594  ±  >16.0937)(5  ±  19.2338) 
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Validation  of  the  LTI  Transfer  Functions 

The  following  figures  validate  the  set  P  of  LTI  plants  used  to  represent  the 
nonlinear  YF-16.  Each  plot  on  the  left  contains  two  outputs.  One  being  the  actual 
output  of  the  YF-16  simulator,  the  other  the  output  of  the  LTI  plant  for  the  same 
input  signal  to  the  elevator’s  actuators.  Each  plot  on  the  right  is  the  error  between 
the  LTI  response  and  that  of  the  YF'16. 
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Figure  A. 6.  Outputs  of  LTI  Plant  ^6  and  of  YF-16  for  the  Same  Input 


Figure  A. 7.  Outputs  of  LTI  Plant  ^7  and  of  YF-16  for  the  Same  Input 
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Figure  A. 12.  Outputs  of  LTI  Plant  ^12  and  of  YF-16  for  the  Same  Input 


Figure  A. 13.  Outputs  of  LTI  Plant  #13  and  of  YF-16  for  the  Same  Input 
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Appendix  B.  Boundary  Data 


This  appendix  contains  the  data  required  to  generate  the  tracking  bounds 
on  the  nonainal  loop  transmission  of  Chapter  IV.  The  additional  boundary  on  the 
2g  responses  is  presented  along  with  the  tabulated  allowable  magnitude  variation 
for  Ig  and  2g  command  bounds.  The  plant  templates  representing  P{j^)  for  the 
frequencies  used  in  generating  the  boundaries  are  also  included. 

Frequency  Domain  Tracking  Bounds: 

It  is  shown  in  Figure  B.l  that  the  2g  responses  used  in  generating  the  equiv 
alent  LTI  plant  set  do  not  fully  populate  the  envelope  of  acceptable  responses.  In 
order  to  make  the  amount  of  uncertainty  more  realistic  an  additional  upper  bound 
is  placed  on  these  outputs.  The  2g  responses  superimposed  upon  the  2g  boundaries 
^  ^  are  shown  in  Figure  B.2,  the  lower  boundary  is  the  same  for  Ig  and  2g  responses. 

The  transfer  functions  representing  the  envelope  bounds  are  displayed  below 
while  the  allowable  magnitude  variations  vs.  frequency  are  presented  in  Tables  B.l 
and  B.2. 

IG  Upper  fit 

15.8965s  +  14.3414  _  15.8965(s  +  0.9022) 


s2  +  8.0691s  +  14.3165 

(s-f  5.4349) (s  +  2.6342) 

2G  Upper  fit 

7.4435s  +  11.3600 

7.4435(s  +  1.5262) 

s2  +  4.2224s  +  11.1403 

(s  + 2.1112  ±2.5852j) 

Lower  fit  for  both  cases 

19.2109 

19.2109 

s2  +  13.6822s  +  19.2703 

(s+  12.088)  (s+  1.5942) 
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Figure  B.l.  2g  Responses  within  Original  Ig  Envelope 


Figure  B.2.  2g  Responses  within  2g  Envelope 
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0.0604 


0.9880 


2.7661 


5.7770 


6.0609 


5.7344 


3.9240 


2.6459 


-1.5395 


-2.3699 


Table  B.l.  Allowed  Variation  in  Magnitude  vs.  Frequency  (ig  Bounds) 


w  (Rad/sec)  Upper  Mag  dB  Lower  Mag  db  |  S/t(u;)  dB 


Lower  Mag  db 

^ii(w)  dB 

-0.0442 

0.1045 

-0.4417 

1.4298 

-1.4975 

4.2636 

-5.5986 

11.3756 

-6.8583 

12.9192 

-11.0619 

16.7962 

-15.7846 

19.7086 

-18.3495 

20.9954 

-23.5948 

23.4470 

-26.1905 

24.6510 

-27.7500 

25.3801 

-42.5389 

32.5239 

3. 

754 

0.1898 


0.6493 


1.8435 


5.4556 


5.8 


-0.2971 


-2.3227 


-5.9656 


-0.0442 


-0.4417 


-1.4975 


-5.5986 


-6.8583 


0.2340 


1.0911 


3.3409 


11.0542 


12.6599 


-15.7846 


-18.3495 


-23.5948 


15.4876 


16.0268 


17.6292 


18.0  i 

-7.5839 

-26.1905 

18.6066 

20.0 

-8.5146  1 

-27.7500 

19.2354 

50.0 

-16.5322 

-42.5389 

26.0068 

Table  B.2.  Allowed  Variation  in  Magnitude  vs.  Frequency  (2g  Bounds) 


Plant  Templates 

The  plant  templates  for  the  frequencies  in  Tables  B.l  and  B.2  are  shown  in 
Figures  B.3  through  B.5  and  are  all  the  same  scale. 
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Vertical  division  is  2  dB 
Horizontal  division  is  10° 


Figure  B.3.  Plant  Templates  for  0.1  <  w  <  5  rad/i 


25  rad/s 


50  rad/sec 


100  rad/sec 


250  rad/sec 


500  rad/sec 


Vertical  division  is  2  dB 
Horizontal  division  is  10° 


Figure  B.5.  Plant  Templates  for  25  <  w  <  500  rad/sec 


Appendix  C.  Simulation  of  the  Stability  Augmentation  System 


This  appendix  presents  the  scheduling  of  the  prefilter  gain  Kf  vs.  com¬ 
manded  input,  along  with  the  normalized  responses  of  the  SAS  at  the  flight  con¬ 
ditions  listed  in  Chapter  V. 

Generation  of  the  Nonlinear  Pre filter 

As  described  in  Chapter  V,  a  nonlinear  prefilter  is  designed  to  better  linearize 
the  final  values  of  the  response  at  the  nominal  point  of  0.9M  and  20,000  ft.  The 
discrete  data  points  for  A/  in  Figure  C.l  are  obtained  by  applying  step  commands 
of  strength  O.lg  to  3.0  g  every  O.lg,  and  at  4g  to  8g  every  Ig.  These  points  are  then 
approximated  by  the  continuous  line  in  the  figure.  This  peace-wise  gain  function 
is  then  programed  in  Fortran  for  use  in  the  nonlinear  YF-16  simulations. 


Figure  C.l.  Scheduling  of  the  Prefilter  Gain  Kf 


{ 


Results  of  Expanding  the  Envelope 

The  following  figures  are  representative  of  the  nonlinear  simulations  run  on 
the  inner  loop  SAS.  The  25  flight  conditions  shown  are  the  same  as  the  tabulated 
data  used  to  present  the  acceptable  envelope  of  Chapter  V.  The  responses  in 
each  figure  are  shown  superimposed  upon  the  normalized  tracking  bounds.  The 
responses  shown  are  only  up  to  the  maximum  command  resulting  in  an  acceptable 
response  in  increments  of  Ig.  When  this  maximum  command  is  not  an  integer, 
the  maximum  acceptable  response  is  also  included.  For  example  Figure  C.2  con¬ 
tains  the  normalized  responses  at  5,000  ft  and  0.5M  with  a  maximum  acceptable 
command  of  7.3g,  the  traces  are  those  for  commands  of  Ig  ,2g  ,3g  ,4g  ,5g  ,6g  ,7g 


Figure  C.2.  Normalized  Response  at  5,000  ft  0.5M 
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TIME 


Figure  C.IO.  Normalized  Response  at  10,000  ft  0.8M 
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Figure  C.14.  Normalized  Response  at  15,000  ft  0.7M 
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Figure  C.21.  Normalized  Response  at  20,000  ft  0.9M 


Figure  C.22.  Normalized  Response  at  30,000  ft  0.5M 
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Figure  C.25.  Normalized  Response  at  30,000  ft  0.8M 
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Response  at  Sf 


Figure  C.26.  Normalized 
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Appendix  D.  Fourth  Order  Fade  Approximation 


The  following  Fade  method  is  taken  directly  from  an  appendix  of  a  digital 
control  systems  text  [11]  and  applied  to  the  function  e’’*  to  obtain  an  approximation 
of  a  pure  time  delay. 


The  General  Fade  Method 

The  Fade  method  produces  rational-fraction  approximations  of  the  form 

Am(l) 


for  the  differentiable  power  series 


m  =  0, 1, 2, . . . 


(D.l) 


d{x)  =  Co  +  CjX  +  cjx^  +  Csx^  +  ...  [D.2) 

where  co  7^  0  and  F{x)  «  g{x).  Am{x)  and  Bm{r)  are  polynomials  of  the  form 

An{x)  =  Co  +  aix  -f-  a2X^  +  H - h  (F.3) 

Bm(x)  =  60  +  biX  -h  -h  631^  -!-•••  +  bmX^  {F.4:) 

where,  for  this  discussion,  the  coefficients  o,  and  bj  are  finite  values  other  than 
zero. 

Let 

r 

Oq 
fll 

a  =  1  aj 


(D.5) 


b  = 


bo 

bi 

62 


(P.6) 
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^p.q  — 


Cg 

Cg~l 

Cg-p 

^9+1 

C9-P+I 

Cj+p-l 

C9+P-2 

Cg-1 

{D.7) 


where  the  elements  of  Eq  (D.l)  whose  subscripts  are  negative,  that  is  i  —  j  <  0, 
are  all  zero. 

The  coefficients  of  Eqs  (D.3)  and  (D.4)  may  be  solved  from  the  following 
matrix  equations: 

Cjn+l.ob  =  a  [D.Q) 

The  matrices  Cni,ni+i  and  Cm+i.o  are  obtained  from  Eq  (D.7),  where  for  the 
former  p  =  m  and  q  =  m  +  1  and  for  the  latter  p  =  m  +  1  and  q  =  0.  Note  that  a 
and  b  are  (m  +  1)  x  1  vectors,  Cni,m+i  is  an  m  x  (m  +  1)  matrix  and  CV+i.o  is  an 
(m  +  l)  X  (m  +  1)  matrix. 


The  Approximation 

The  power  series  of  e*  where  i  =  rs  is  expanded  in  a  MacLaurin  series  about 
zero,  resulting  in 

d(i)  =  e*  =  l  +  i+^  +  |j-  +  ^  +  ---  (D.IO) 

d{x)  is  to  be  approximated  by  the  fourth  order  rational  fraction 

.  _  Op  -f  Oil  +  a2X^  +  031^  +  ajX*  f  D 

6o  +  bix  +  b2x'^  +  i>3i^  +  biX* 

where  m  =  1,  cp  =  ci  =  1,  C2  =  j,C3  =  |,  etc.  From  Eq  (D.7)  the  following  are 
obtained: 


C4,S  = 


C5 

C4 

C3 

C2 

Cl 

C6 

C5 

C4 

C3 

C2 

c? 

ce 

Cs 

C4 

C3 

C8 

Cl 

Ce 

C5 

C4 

D-2 


Cs.o  = 


Co 

0 

0 

0 

0 

Cl 

Co 

0 

0 

0 

C2 

Cl 

Co 

0 

0 

C3 

C2 

Cl 

Co 

0 

C4 

C3 

C2 

Cl 

Co 

Inserting  these  matrices  into  Eqs  (D.8)  and  (D.9)  results  in  for  Eq  (D.8) 


= 


1 

1 

1 

1 

1 

S! 

24 

6 

2 

1 

1 

1 

1 

1 

6! 

S! 

24 

6 

2 

1 

1 

1 

1 

1 

7! 

6! 

5! 

24 

6 

1 

1 

1 

1 

1 

.  8! 

7! 

6! 

5! 

24 

bo 

bi 

62 

bz 

64 


=  0 


There  are  only  4  equations  and  5  unknowns,  setting  i>o  =  1  Eq  (D.12)  becomes 


5  20  60  120 

6  30  120  360 

7  42  210  840 

8  56  336  1680 


bo 

-1 

61 

-1 

62 

II 

-1 

^3 

-1 

——  -J 

-1 

(D.13) 


Solving  for  b  gives 


_3 

28 

84 

1 

1680 
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Inserting  the  values  of  Cm+i.o  and  b  into  Eq  (D.9)  results  in 


1 

0 

0 

0 

0 

1 

Oo 

1 

1 

0 

0 

0 

-1 

2 

ai 

1 

2 

1 

1 

0 

0 

3 

28 

= 

dj 

1 

6 

1 

2 

1 

1 

0 

-1 

84 

03 

1 

24 

1 

6 

1 

2 

1 

1 

1 

1680 

04 

Solving  for  a  gives 


a  = 


1 

2 

28 

J_ 

84 


1680 


Inserting  the  values  of  a,  b  and  x  —  rs  into  Eq  (D.ll)  yields 


F(x)  = 


2953 


.  840  c  u- 

I  7^0  i  ^4 


54  _  ^^3  +  ^^2  -ms  + 


1680 


{P.14) 


For  the  pilot  model  of  Chapter  VI  the  time  delay  of  r  =  -0.3  gives  the  approxi¬ 
mation  used  in  the  nonlinear  simulations  as 


s*  -  66.66752  +  20005*  -  3.111  X  10^5  +  2.0741  x  10^ 
s*  +  66.66753  ^  2OOO52  +  3.111  X  10<5  +  2.0741  X  105 


(5  -  14.03  ±  il7.7)(5  -  19.31  ±  j5.78) 
(5  +  14.03  ±  j17.7)(5  +  19.31  ±  j5.78) 


(D.15) 
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